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a.  I 

ABSTRACT 


The  Lorentz  force  equation  governing  the  interaction  between  radiation  and 
charged  particles  in  the  presence  of  a  static  magnetic  field  can  be  reduced  to  the  form 
of  the  simple  pendulum  equation,  when  applied  to  the  Free  Electron  Laser  (FEL)  and  to 
the  amplification  of  whistler  waves  in  the  magnetosphere. 

The  first  topic  in  this  thesis  is  the  start-up  of  an  FEL  with  a  compact  design  having 
a  small,  spatial  undulator  period,  many  periods,  and  a  small  undulator  field  strength. 
Under  these  conditions,  the  number  of  photons  spontaneously  emitted  into  the 
resonator  mode  each  pass  is  so  small  that  a  classical  field  cannot  be  established. 
Quantum  fluctuations  affect  the  wave-particle  interaction  and  therefore  diminish  weak- 
field  gain.  An  FEL  start-up  condition,  several  orders  of  magnitude  more  restrictive  than 
determined  by  quantum  recoil  alone,  is  derived  from  the  photon  rate  equation. 

The  next  topic  is  the  evolution  of  a  weak,  classical  radiation  field  driven  by  a  long 
electron  pulse.  Assuming  low  gain,  it  is  possible  to  reduce  the  coupled  optical  wave 
and  electron  pendulum  equation  to  a  simple,  first-order,  partial  differential  equation 
including  electron  dynamics,  slippage,  desynchronism,  and  resonator  loss.  This 
analytical  method  is  used  to  study  coherence  development,  and  the  effects  of  varying 
the  electron  beam  pulse  shape. 

Next,  the  Boeing  Average  Power  Laser  Experiment  (APLE)  is  investigated  using 
two  and  three  dimensional  numerical  simulations.  Various  methods  of  reducing  the 
electron-radiation  coupling  in  the  oscillator  are  examined  in  order  to  minimize  the 
electron  beam  energy  spread  introduced  into  the  FEL  amplifier.  An  alternative  for 
achieving  high  power  makes  use  of  the  oscillator  alone,  but  with  a  tapered  undulator  to 
increase  the  electron-radiation  coupling. 

The  final  topic  is  charged  particles,  following  helical  trajectories,  traveling  in  the 
magnetosphere  between  mirror  points  along  the  earth's  magnetic  field  lines.  When 
electrons  interact  with  a  whistler  wave  near  the  geomagnetic  equator,  their  motion  can 
again  be  described  by  the  simple  pendulum  equation.  The  inhomogeneity  of  the 
earth's  magnetic  field  near  the  geomagnetic  equator  is  analogous  to  the  field  in  a 
tapered  undulator.  Numerical  phase-space  simulations  are  developed  to  investigate 
the  interaction. 


IV 


ducLey  KNtrnte1 

NAvkL  POsfJMyRffiPtfghHnn, 


REPORT  DOCUMENTATION  PAGE 


:  reporting  burden  for  this  collection  of  information  is  estimated  to  average  I  hour  per  response,  including  the  t'ftBilWM>WBPlPCr''Jcy*nA  se^£WHljf'!ls,.Cl<h^  sourceV 
ring  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information    Send  commenY*Vwar«rl^n^RjrdeW4^maH<o^Ti(J(TMrlaV)Jct  of  this 


Public  i 

gatheni  w 

collection  of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services.  Directorate  tor  information  Operations  and  Reports.  1215  JeHerson 

Davis  Highway.  Suite  1204.  Arlington.  VA  22202-4302  and  to  the  Off  ice  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188).  Washington.  DC  20S03 


1.   AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

June    19  93 


3.  REPORT  TYPE    AND  DATES  COVERED 

Dissertation 


4.  TITLE  AND  SUBTITLE 

A  Study  of  the  Amplification  of  Laser  and  VLF 
Waves  Using  a  Simple  Pendulum  Model 


6.  AUTHOR(S) 


Keith  Alan  Sturgess 


5.    FUNDING  NUMBERS 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 
Department  of  Physics 
Monterey,  California  93943 


PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


Approved  for  Public  Release;  distribution  is  unlimited. 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words)    m,  T  ^  ., 

The  Lorentz  force  equation  governing  the  interactio 
between  radiation  and  charged  particles  in  the  presence  of  a  static 
magnetic  field  can  be  reduced  to  the  form  of  the  simple  pendulum 
equation,  when  applied  to  the  Free  Electron  Laser  (FEL)  and  to  the 
amplification  of  whistler  waves  in  the  magnetosphere. 

The  first  topic  in  this  thesis  is  the  start-up  of  an  FEL  with  a 
compact  design  having  a  small,  spatial  undulator  period,  many  periods, 
and  a  small  undulator  field  strength.   Under  these  conditions,  the 
number  of  photons  spontaneously  emitted  into  the  resonator  mode  each 
pass  is  so  small  that  a  classical  field  cannot  be  established.   Quantum 
fluctuations  affect  the  wave-particle  interaction  and  therefore  diminis 
weak  field  gain.   An  FEL  start-up  condition,  several  orders  of  magnitud 
more  restrictive  than  determined  by  quantum  recoil  alone,  is  derived 
from  the  photon  rate  equation. 

The  next  topic  is  the  evolution  of  a  weak,  classical  radiation 
field  driven  by  a  long  electron  pulse.   Assuming  low  gain,  it  is 


14.  SUBJECT  TERMS 


Free  Electron  Laser,  FEL,  undulator,  start-up, 
quantum  fluctuations,  longitudinal  modes,  whistler 

wav^s.     ppnrjnliim    eqnaHnnr     niimpriral     simii1at; 


fl 


i  nn 


15.  NUMBER  OF  PAGES 

103 


16.  PRICE  CODE 


17.    SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.    SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF  ABSTRACT 

UL 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev    2-89) 


REPORT  DOCUMENTATION  PAGE,  SF  2  98 


13.   Continued 


possible  to  reduce  the  coupled  optical  wave  and  electron 
pendulum  equation  to  a  simple,  first-order,  partial 
differential  equation  including  electron  dynamics,  slippage, 
desynchronism,  and  resonator  losses.   This  analytical 
method  is  used  to  study  coherence  development,  and  the 
effects  of  varying  the  electron  beam  pulse  shape. 

Next,  the  Boeing  Average  Power  Laser  Experiment  (APLE) 
is  investigated  using  two  and  three  dimensional  numerical 
simulations.   Various  methods  of  reducing  the  electron- 
radiation  coupling  in  the  oscillator  are  examined  in  order 
to  minimize  the  electron  beam  energy  spread  introduced 
into  the  FEL  amplifier.   An  alternative  for  achieving  high 
power  makes  use  of  the  oscillator  alone,  but  with  a  tapered 
undulator  to  increase  the  electron-radiation  coupling. 

The  final  topic  is  charges  particles,  following  helical 
trajectories,  traveling  in  the  magnetosphere  between  mirror 
points  along  the  earth's  magnetic  field  lines.   When  electrons 
interact  with  a  whistler  wave  near  the  geomagnetic  equator, 
their  motion  can  again  be  described  by  the  simple  pendulum 
equation.   The  inhomogeneity  of  the  earth's  magnetic  field 
near  the  geomagnetic  equator  is  analogous  to  the  field  in  a 
tapered  undulator.   Numerical  phase-space  simulations  are 
developed  to  investigate  the  interaction. 


11 


Table  of  Contents 


I.  INTRODUCTION 1 

II.  FREE  ELECTRON  LASER  THEORY 6 

A.  INTRODUCTION 6 

B.  ELECTRON  DYNAMICS 9 

C.  THE  SELF-CONSISTENT  WAVE  EQUATION 12 

D.  FEL  DIMENSIONLESS  PARAMETERS 16 

III.  QUANTUM  LIMITATION  ON  FREE  ELECTRON  LASER  START-UP 18 

A.  INTRODUCTION 18 

B.  COMPACT  FEL  DESCRIPTION 19 

C.  START-UP  OF  THE  COMPACT  FEL 20 

D.  ANALYTICAL  RESULTS 26 

E.  CONCLUSIONS  AND  FURTHER  RESEARCH  SUGGESTIONS 27 

IV.  ANALYTICAL  DESCRIPTION  OF  LONGITUDINAL  MODES 29 

A.  INTRODUCTION 29 

B.  MATHEMATICAL  DEVELOPMENT 30 

C.  ELECTRON  PULSE  EFFECTS 34 

1.  Continuous  Electron  Beam 36 

2.  Parabolic  Electron  Beam 38 

3.  Triangular  Electron  Pulse 44 

D.  COHERENCE  DEVELOPMENT 48 

E.  CONCLUSIONS  AND  FURTHER  RESEARCH  SUGGESTIONS 51 


V 


V.  SHORT  PULSE  EVOLUTION  IN  THE  BOEING  APLE  OSCILLATOR 53 

A.  INTRODUCTION 53 

B.  OSCILLATOR  DESCRIPTION 54 

C.  DIFFRACTION  EFFECTS 56 

D.  PULSE  EFFECTS 57 

1.  No  Energy  Slew  in  the  Undulator 57 

2.  Energy  Slew  in  the  Oscillator „ 60 

3.  High  Power  Oscillator  Experiment 62 

E.  CONCLUSIONS  AND  FURTHER  RESEARCH  SUGGESTIONS 64 

VI.  WHISTLER  WAVE  AMPLIFICATION  IN  THE  MAGNETOSPHERE 66 

A.  INTRODUCTION 66 

1.  Description  of  Whistler  Waves 66 

2.  Description  of  the  Magnetosphere.. 67 

B.  WAVE  PARTICLE  INTERACTIONS 69 

1.  The  Cold  Background  Plasma 69 

2.  The  Whistler  Wave  Equation 72 

3.  The  Hot  Electron  Equations  of  Motion 75 

C.  NUMERICAL  SIMULATIONS 81 

D.  CONCLUSIONS  AND  FURTHER  RESEARCH  SUGGESTIONS  .....  88 

LIST  OF  REFERENCES 90 

INITIAL  DISTRIBUTION  LIST...... 93 


VI 


ACKNOWLEDGEMENT 


As  with  any  significant  effort,  there  are  many  people  directly  and  indirectly 
responsible  for  the  success  of  this  dissertation.  I  would  like  to  first  thank  Bill 
Colson,  who  took  me  in,  gave  me  interesting  projects  to  work  on,  and  stood 
behind  me  and  pushed  when  the  going  got  tough  down  the  home  stretch. 
Second,  my  thanks  go  out  to  Professors  Armstead  and  Neighbours  who  gave 
valuable  input  into  the  writing  of  this  dissertation,  and  whose  assistance  in 
preparing  for  the  oral  exams  was  essential. 

There  is  no  way  to  quantify  all  of  the  help  my  wife,  Lynn,  has  been  to  me 
during  the  writing  of  this  dissertation.  She  was  my  second  reader  for  spelling 
and  grammar.  She  made  all  of  the  changes  for  a  rewrite  while  I  was  away  on 
TDY  in  Kansas  and  Virginia.  She  accompanied  me  on  several  research  and 
conference  trips  and  probably  knows  more  about  free  electron  lasers  than  some 
physicists.  Her  love  and  understanding  were  a  constant  source  of  unwavering 
support  that  helped  me  get  through  some  difficult  days  and  weeks.  I  only  hope  I 
support  her  as  well  as  she  did  me  when  she  is  working  toward  her  doctorate. 


vxx 


I.  INTRODUCTION 


A  free  electron  laser  (FEL)  oscillator  uses  a  beam  of  relativistic  electrons 
passing  through  a  periodic  undulator  to  produce  and  amplify  light  stored  within  a 
resonator.  Since  this  light  is  a  product  of  the  electron  -  undulator  interaction,  its 
properties  (i.e.  wavelength  and  power)  can  be  changed  by  altering  the  properties 
of  the  electron  beam  and  the  undulator.  The  ease  of  tunability  of  the  FEL  over  a 
wide  range  of  wavelengths  is  beginning  to  make  the  FEL  a  sought  after  piece  of 
equipment  by  many  disciplines  in  universities,  research  labs,  and  industry.  The 
promise  of  high  efficiency  and  high  average  power  has  made  the  FEL  a 
research  priority  of  the  armed  forces  as  they  seek  a  speed-of-light  weapon. 

The  basic  physics  of  an  FEL  is  presented  in  Chapter  II  as  a  starting  point  for 
the  research  chapters  that  follow.  The  FEL  theory  presented  here  uses 
Maxwells'  field  equations  to  describe  the  evolution  of  the  light  in  response  to  the 
bunching  electrons,  and  self-consistently  uses  the  pendulum  equation  to 
describe  the  evolution  of  the  bunching  electrons  in  response  to  the  light  [1]. 

Chapter  III  examines  FEL  physics  during  the  first  several  passes  of  the 
electrons  through  a  compact  undulator  that  uses  many  short  periods  and  a  small 
undulator  parameter.  Spontaneous  emission  is  introduced,  and  the  effects  of  a 
small  number  of  photons  in  the  resonator  causes  statistical  fluctuations  of  the 
photon  number.  These  quantum  fluctuations,  if  large  enough,  can  significantly 
degrade  the  classical,  small-signal  gain  of  the  FEL  causing  slow  start-up  or 
possibly  the  failure  of  the  FEL  to  start. 

New  and  original  research  concerning  the  problem  of  quantum  limitations  on 
the  start-up  of  FELs  is  presented  in  Chapter  III.  The  previously  derived  equation 
for  gain  including  photon  fluctuations  is  incorporated  into  the  photon  rate 
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equation  to  obtain  a  new  start-up  condition.  The  photon  rate  equation  is  then 
solved  analytically  using  the  new  expression  for  gain.  This  and  other  analytical 
results  are  compared  with  those  obtained  by  other  authors  using  a  completely 
quantum  mechanical  method.  Also,  the  new  expression  for  gain  is  incorporated 
into  a  numerical  simulation  to  solve  for  optical  mode  evolution. 

Chapter  IV  follows  the  evolution  of  long  optical  pulses  using  a  partial- 
differential  equation  which  includes  effects  of  desynchronism  and  resonator 
losses.  Desynchronism  describes  the  misalignment  between  the  optical  and 
electron  pulses  on  successive  passes  through  the  resonator.  Analytical 
expressions  are  derived  for  the  optical  wave  evolution  assuming  a  long  electron 
pulse,  while  the  FEL  operates  in  the  small  signal,  low  gain  regime.  The  first- 
order,  partial-differential  equation  is  then  used  to  examine  the  effects  of  the 
electron  pulse  shape  on  the  final  optical  pulse.  Also,  the  coherence 
development  from  a  random  incoherent  optical  pulse  to  a  final  coherent  optical 
pulse  is  examined. 

A  partial-differential  equation  describing  optical  pulse  evolution  is  used  for 
new  and  original  research  in  Chapter  IV.  The  single-pass,  first-order  partial- 
differential  equation  is  extended  from  one  to  multiple  passes  including 
desynchronism  and  resonator  losses.  The  partial-differential  equation  is  solved 
analytically  for  the  case  of  a  continuous  electron  beam,  a  parabolic  electron 
pulse  shape,  and  a  Gaussian  electron  pulse.  A  desynchronism  curve  showing 
the  variation  of  optical  power  with  desynchronism  is  then  created  from  the 
analytical  solution  for  a  parabolic  electron  pulse  shape  and  compared  to 
desynchronism  curves  produced  numerically,  experimentally,  and  with  analytical 
methods  used  by  other  authors.  It  is  shown  that  exact  synchronism  gives  zero 
optical  power  in  steady-state.  The  case  of  a  triangular  electron  pulse  is  solved 
analytically.    A    numerical  FEL  simulation  is  modified  to  solve  the  single-pass 


partial-differential  equation  and  used  to  study  coherence  evolution.  Finally,  a 
numerical  code  is  developed  to  compare  the  analytical  solution  and  the 
numerical  FEL  simulation  for  the  cases  of  a  parabolic,  triangular  or  continuous 
electron  pulse. 

Only  weak  optical  fields  are  studied  through  Chapter  IV.  Chapter  V 
examines  the  Boeing/LANL  high  Average  Power  Laser  Experiment  (APLE). 
Since  the  proposed  design  of  APLE  has  a  single  electron  beam  powering  both 
an  oscillator  and  amplifier  FEL  in  series,  it  is  important  that  the  oscillator  FEL  not 
induce  a  large  energy  spread  in  the  electron  beam  prior  to  the  amplifier  FEL. 
Chapter  V  examines  methods  of  reducing  the  electron  -  light  coupling  in  the 
oscillator.  One  method  of  reducing  the  coupling  calls  for  inducing  a  ramp,  or  a 
slew,  in  energy  along  the  electron  pulse.  Running  the  oscillator  alone  as  a  high 
average  power  FEL  is  also  examined. 

The  investigation  into  methods  of  reducing  the  electron  beam  energy  spread 
induced  by  the  FEL  oscillator  at  the  Boeing  APLE  project  is  new  and  original 
research.  A  numerical  FEL  simulation  is  extended  to  accurately  calculate  the 
optical  field  evolution  in  the  presence  of  large  desynchronism.  Using  this 
modification,  a  desynchronism  curve  is  presented  for  the  APLE  oscillator.  The 
desynchronism  that  best  reduces  the  electron  energy  spread  for  the  APLE 
oscillator  is  found.  The  simulation  is  modified  further  to  include  the  effects  of 
energy  and  timing  jitter  of  the  electron  pulse  from  the  accelerator.  These  effects 
are  explored  for  large  and  small  values  of  desynchronism.  The  simulation  is 
modified  again  to  include  an  energy  slew  across  the  electron  micropulse  in  order 
to  show  that  the  energy  spread  induced  in  the  electron  beam  is  more  than 
Boeing  researchers  had  expected.  The  final  APLE  topic  explores  the  tapered 
oscillator  design  as  a  high  average  power  device.  Conditions  are  obtained  for 
achieving  the  optimum  power  from  the  FEL  oscillator  alone  with  no  amplifier. 


Finally,  Chapter  VI  examines  whistler  wave  amplification  employing  the 
similar  theoretical  and  simulation  methods  used  to  understand  the  FEL.  The 
whistler  and  FEL  both  have  electrons  in  helical  trajectories  coupling  to  and 
amplifying  a  transverse  electromagnetic  wave.  The  wave  and  electron 
equations  of  motion  are  developed  for  the  whistler  interaction  including  the 
inhomogeneity  due  to  the  variation  in  magnetic  field  strength  along  a  field  line  of 
the  earth's  magnetic  field.  Simulations  are  then  developed  from  these  equations 
in  order  to  describe  stimulated  VLF  emission. 

Chapter  IV  includes  new  and  original  research  on  whistler  wave 
amplification  using  a  wave  equation  driven  by  the  pendulum  equation.  As  a 
member  of  a  team,  I  participated  in  the  derivation  of  the  homogeneous  and 
inhomogeneous  wave  equation  and  electron  equations  of  motion.  Next,  the 
team  developed  numerical  simulations  of  VLF  wave  amplification.  Using  the 
numerical  simulations,  the  effects  of  homogeneous  and  inhomogeneous 
magnetic  fields  are  explored,  as  well  as  the  effects  of  a  spread  in  phase  velocity. 
The  phase-  space  evolution  of  the  electrons  and  the  growth  of  the  optical  pulse 
power  and  phase  are  displayed  graphically  in  order  to  better  understand  their 
coupled  dynamics. 

Much  of  the  research  in  this  thesis  has  been  presented  at  scientific 
conferences,  published,  or  submitted  for  publication.  The  research  in  Chapter  III 
will  be  submitted  to  Physical  Review  Letters  under  the  title  "A  Quantum 
Limitation  on  the  Start-up  of  Compact  Free  Electron  Lasers"  by  K.  A.  Sturgess 
and  W.  B.  Colson.  The  research  in  Chapter  IV  was  presented  at  the  Thirteenth 
International  Free  Electron  Laser  Conference  in  Sante  Fe,  New  Mexico  (August 
1991),  and  resulted  in  the  publication  of  a  paper  titled  "An  Analytical  Solution  to 
Longitudinal  Modes  in  Free  Electron  Lasers"  by  K.  A.  Sturgess  and  W.  B. 
Colson  in  Nuclear  Instruments  and  Methods  in  Physics  Research,  A318,  pages 


576-581,  (1992).  The  research  in  Chapter  V  was  presented  at  the  Fourteenth 
International  Free  Electron  Laser  Conference  in  Kobe,  Japan,  (August  1992) 
and  accepted  for  publication  in  Nuclear  Instruments  and  Methods  in  Physics 
Research  with  the  title  "Short  Pulse  Evolution  in  the  Boeing/LANL  APLE 
Oscillator"  by  K.  Sturgess,  D.  J.  Frost  and  W.  B.  Colson.  The  research  in 
Chapter  VI  was  presented  at  the  Fourteenth  International  Free  Electron  Laser 
Conference  in  Kobe,  Japan,  (August  1992),  and  accepted  for  publication  in 
Nuclear  Instruments  and  Methods  in  Physics  Research  with  the  title  "Analogies 
Between  the  Free  Electron  Laser  and  Whistler  Wave  Amplification  in  the 
Magnetosphere"  by  W.  B.  Colson,  B.  Gately,  D.  L.  Caudle,  and  K.  Sturgess. 


I.  FREE  ELECTRON  LASER  THEORY 


A.  INTRODUCTION 

An  operating  Free  Electron  Laser  (FEL)  consists  of  an  accelerator  capable 
of  driving  electrons  to  relativistic  energies,  an  undulator  composed  of  a 
transverse,  periodic  magnetic  field  that  "wiggles"  the  electrons,  and  optical 
devices  necessary  to  collect  and  transport  the  amplified  radiation.  The 
periodically  deflected  relativistic  electrons  radiate  spontaneously,  and  in 
addition,  couple  to  the  radiation  store  in  an  optical  resonator  leading  to  coherent 
stimulated  emission.  The  idea  of  an  FEL  was  conceptually  introduced  in  1971 
by  J.M.J  Madey  [2].  Related  work  which  contributed  to  the  FEL  include  the 
traveling  wave  tube  of  Motz  [3],  and  the  ubitron  developed  by  Philips  [4]. 
Successful  experimental  operation  of  the  FEL  amplifier  in  1976  [5]  and  an  FEL 
oscillator  in  1977  [6]  demonstrated  the  feasibility  and  potential  of  this  new  laser 
device. 

As  mentioned  above,  the  undulator's  role  in  an  FEL  is  to  give  the  relativistic 
electrons  periodic  transverse  deflections  as  they  travel  with  the  light  pulse  down 
the  length  of  the  undulator.  This  deflection,  provided  by  either  permanent 
magnets  or  electromagnets,  allows  the  electrons  to  couple  to,  and  transfer 
energy  with  the  transverse  optical  pulse  along  the  undulator  length,  L  =  N/L0 
where  A0  is  the  undulator  wavelength  and  N  is  the  number  of  periods  in  the 
undulator.  An  important  dimensionless  quantity  characterizing  the  undulator  is 
the  undulator  parameter  K  =  eBA0/2nmc2 ,  where  B  is  the  rms  field  strength 
over  one  period  of  the  undulator,  e  is  the  electron  charge  magnitude,  m  is  the 
mass  of  the  electron,  and  c  is  the  speed  of  light.  In  most  FEL's,  K*0.1->5. 


There  are  two  main  types  of  undulator,  based  on  the  paths  electrons  take 
within  them:  helical  and  linearly  polarized.  A  helical  undulator  produces  helical 
electron  paths,  which  produces  and  couples  to  circularly-polarized  light.  The 
ideal  helical  magnetic  field  is  constant  in  magnitude  along  the  longitudinal  axis 
of  the  undulator,  and  in  rectangular  coordinates  is  given  by  [1] 

BH  =  B[cos{k0z),sin(k0z),0].  (2-1) 

where  B  is  the  peak  magnetic  field  strength  on  this  axis,  and  k0  =  2n/A0  is  the 
undulator  wave  number.   The  helical  undulator  design  focuses  electrons  toward 
the  undulator  axis,    since  the  average  magnetic  field  can  be  expanded  and 
shown  to  increase  off-axis  as 


B  =  B 


k2r2 


\ 


1  +  n°,4+...  ,  (2-2) 


where  r  is  the  distal  distance  from  the  axis  of  the  undulator. 

A  linearly-polarized  undulator  produces  periodic  planar  motion  in  the 
electron  beam  which  in  turn  produces  and  couples  to  linearly-polarized  light. 
The  linear  magnetic  field  has  an  rms  field  strength  of  S/V2  which  is  less  than 
the  helical  case,  and  so,  couples  less  well  to  the  light  and  produces  half  the 
gain.   In  rectangular  coordinates,  the  linear  field  is  given  by 

BL  =  B[0,sin(k0z),0]  (2-3) 

and  provides  magnetic  focusing  in  one  direction  only. 

The  accelerators  used  in  various  FEL  projects  range  from  low  energy  Van 
De  Graff  electrostatic  accelerators  to  high  energy  storage  rings,  including 
intermediate  energy  microtrons  and  linear  accelerators.    They  span  an  energy 


range   giving    Lorentz  factors    7«10-»103,   where    y  1  =  ^1  -  v2/c2  ,      and 


electron  beam  currents  ranging  from  I«1->103A.  The  typical  RF  linac 
produces  electron  micropulses  of  length,  /e  «  mm  -»  cm  ,  separated  by  about 
one  to  several  meters  grouped  together  in  a  long  macropulse  of  1/iS  to  1ms. 
The  transverse  radius  of  the  electron  beam  is  on  the  order  of  millimeters 
resulting  in  an  electron  density  of  p«  1012  ->  1013  cm"3. 

The  electron  beam  quality  is  characterized  by  its  energy  spread  and 
emittance.  Energy  spread  describes  the  spread  of  electron  velocities  in  the 
beam  and  is  characterized  by  Ay/y.  Emittance  is  defined  as  £•«  r0^,  where  r0 
is  the  rms  initial  position  spread  of  the  electrons  in  the  beam,  and  ~60  is  the  rms 
initial  angular  spread  in  the  beam.  Although  the  angular  and  position  spread  of 
the  beam  can  be  varied,  the  emittance  of  the  beam  tends  to  be  a  constant  for  an 
accelerator  at  a  given  energy.  Normalized  emittance,  en  «  ye,  tends  to  be 
constant  along  the  accelerator  while  the  electrons  are  increasing  their  energies. 

When  the  electron  beam  enters  the  undulator,  it  is  desirable  to  choose  the 
angular  and  position  spreads  in  the  beam  so  that  it  does  not  excessively  focus 
or  expand  as  it  travels  down  the  length  of  the  undulator.  For  a  fixed  value  of 
emittance,  a  small  beam  radius  at  one  point  along  the  undulator  results  in  a 
large  angular  spread,  and  consequently  a  large  beam  radius  farther  along  the 
undulator.  A  small  angular  spread  results  in  a  large  beam  radius  everywhere. 
The  condition  that  minimizes  the  beam  radius  along  the  undulator  is  called  the 
'matching  condition,'  and  gives  Kk0rQ  «  y\  [1],  This  condition  determines  an 
electron  beam  radius  of  re  =  >Jyenz/y[27rK  where  A  =  2n/k  is  the  optical 
wavelength.  A  beam  matched  in  this  way  gives  a  phase  velocity  spread  of 
Av=27rf\IKk£/y.  Since  the  phase  velocity  spread  must  be  less  than  the  gain 
band-width,  Av<2;r,  an  upper  limit  on  the  allowable  emittance  is 
(Tax  <yA/2xNK. 
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There  have  been  many  approaches  to  understanding  how  the  FEL  works. 
Among  the  methods  used  to  develop  FEL  theory  has  been  a  quantum 
mechanical  approach,  the  use  of  quantum  electrodynamics  ,  the  application  of 
plasma  physics  theories,  and  the  use  of  plasma  physics  many-particle  computer 
codes  [1].  Early  on  in  the  theoretical  development,  it  was  shown  that  a  classical 
approach  provided  an  accurate  description  of  the  FEL  [7],  The  classical 
approach  using  a  charged-particle  current  driving  Maxwell's  equations  will  be 
presented  here. 


B.  ELECTRON  DYNAMICS 

The  only  significant  force  acting  on  an  electron  as  it  travels  the  length  of  the 
undulator  is  the  Lorentz  force  given  by  the  relativistic  equations  [1] 


d[ip)       -e(£f+^(Br+B„))  P-») 


dt         mc 


%  =  —{p-E,)  (2-5) 

dt       mc 

y-2  =  \-~p~p  (2-6) 

where  /?  =  v/c  is  the  electron's  normalized  velocity,  Bu  =  B[cos(k0z) ,sin{k0z) ,0] 
is  the  helical  undulator  magnetic  field,  Er  =  E[cos  y/,-sin  y/,0]  and 
Br  =  E[sin  y,cos  y/,0]  are  the  optical  electric  and  magnetic  fields, 
y/  -  kz  -  dt  +  0,  k  =  2x/A  is  the  optical  wavenumber,  co  =  ck  is  the  optical 
frequency,  and  <p  is  the  phase  of  the  optical  wave.  Equations  of  motion  for  the 
electrons  are  obtained  by  inserting  the  electric  and  magnetic  fields  into  the 
equations  above.    The  plane  wave  form  for  the  electric  and  magnetic  fields  is 


accurate  near  the  axis  of  the  undulator.  Since  there  are  a  total  of  five  equations 
above,  and  only  four  unknowns,  [r(f),y(f)]f  one  of  the  equations,  the  z 
component  of  equation  (2-4),  will  be  ignored.  Expanding  the  remaining 
equations  yields 

— y~  =  —  \E^-  0z)[cos  v,-sin  v,0]  +  /3zB[-sin(k0z),cos(k0z),0]}  (2-7) 
at  mcl  L  JJ 


dy      -eE 


{pxcosY-pysiny,)  (2-8) 


dt       mc 

r~2  =  i-~0i-%  (2-9) 

where  J}1  =  [/?,,/?y)0]  is  the  transverse  electron  velocity.    Since  the  electrons  in 

an  FEL  are  relativistic,  vz  «c,  we  have   (1-/?J£  «/3zB,  so  that  the  transverse 

optical  force  is  much  smaller  than  the  undulator  field  force.    This  removes  the 

first  term  on  the  right  side  of  equation  (2-7).  The  remaining  term  can  be  directly 

integrated  yielding 

\  =  —  [cos(/c0z),s/h(/c0z),0]  (2-10) 

Y  ■ 

The  constants  of  integration  are  taken  to  be  zero  representing  an  electron  beam 

that  is  injected  perfectly  into  helical  orbits. 

The  transverse  amplitude  of  the  electrons  can  be  found  from  integrating 

equation  (2-10)  further.    First,  assume  that  y(t)  is  a  constant,  which  is  a  good 

approximation  for  low  efficiency  FELs.    Next,  assume  that  z  &fi0ct,  so  that  the 

electron  oscillation  frequency  in  the  undulator  is  co^=kJ3Qc.    The  arguments 

within  the  trigonometric  functions  of  Equation  (2-10)  can  then  be  rewritten  as 

k0z  »co0t .   Integrating  now  leads  to 

x,  * 


±  *^[-sin(wQt),cos(co0t),0]  (2-11) 

Z7ryL  J 
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which  yields  a  transverse  amplitude  of  oscillation,  Ax±  «  KA0/27ry. 

Continuing  with  pendulum  equation  derivation,  substitute  equation  (2-10) 
into  the  energy  equation  (2-8).  Using  a  trigonometric  identity  leads  to 

^  =  —  cos((/c0  +  k)z  -cot  +</>).  (2-12) 

at       ymc 

The  electron's  phase  is  identified  in  equation  (2-12)  as  £=(k  +  k0)z-cot.  We 
see  from  equation  (2-12)  that  electron  phases  such  that  -nj2<  £+</><  /r/2  will 
gain  energy  from  the  radiation,  while  electron  phases  such  that 
n/2  <  £+  <f>  <  Zk  2  will  lose  energy  to  the  radiation. 

Differentiating  the  electron  phase  velocity  twice  yields  £=  (k  +  k0)/32c. 
Inserting  the  magnitude  of  equation  (2-10)  into  the  definition  of  y  from  equation 
(2-9)  and  differentiating  with  respect  to  time  yields  y/y  =  y2fljlz .  Combining 
these  equations  with  equation  (2-12)  gives  the  pendulum  equation 

■i=2el^KEcos(i+^  {2_13) 

y  m 

which  governs  the  electron's  motion  in  the  presence  of  the  combined  radiation 
and  undulator  fields.  Defining  the  dimensionless  time  as  z=ct/L  (so  that 
r  =  0  ^>1  as  the  electrons  travel  from  the  beginning  to  the  end  of  the  undulator) 
normalizes  equation  (2-13)  as 

°^r)  =  ^2e^KEcos^+    j  =  |a|cos(c+  ^  (2.14) 

where  (..)  =  d2{.)t  dr2  =  {Lc)  d2(..)/dt2  and  |a|  =  4xNeEKL/ylmc2  is  the 
dimensionless  optical  field  strength.  The  relativistic  Lorentz  factor  y  in  equation 
(2-13)  has  been  replaced  by  its  initial  value  y0  in  equation  (2-14)  by  assuming 
low  efficiency:   the  energy  changes  of  the  electrons  are  insignificant  during  the 
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interaction  time  r.  The  electron  pendulum  equation  for  a  helical  undulator  (2- 
14)  is  appropriate  for  low  and  high  gain  and  correctly  follows  the  development  of 
electron  bunching.  For  the  linear  undulator,  \a\  should  be  modified  by  letting 
K->  K(j0(g)-  »/,(£)),  where  J0  and  J,  are  Bessel  functions  and 
£=K2/2(l  +  K2). 

Using  the  expressions  above,   a  dimensionless  phase  velocity  can   be 
defined  by 

v{r)  =  ^=L[(/c  +  /c0  )&-*].  (2-15) 

The  electron  experiences  resonant  undulator  and  radiation  field  forces  when 
v=0.  The  resonant  relationship  between  the  radiation  and  undulator 
wavelengths  is 

^  =  ^0-^-2-  (2"16) 

2/ 

Equation  (2-16)  shows  that  the  operating  wavelength  of  an  FEL  can  be  altered 
by  changing  the  electron  beam  energy  y,  the  undulator  wavelength  A0,  or  the 
undulator  parameter  K. 


C.  THE  SELF-CONSISTENT  WAVE  EQUATION 

While  the  Lorentz  force  law  governs  the  evolution  of  the  electron's  trajectory, 
Maxwell's  equations  govern  the  evolution  of  the  optical  field  in  the  FEL. 
Spontaneous  emission  in  an  FEL  oscillator  grows  to  form  a  classical  field  with  a 
bandwidth  proportional  to  the  inverse  of  the  number  of  periods  in  the  undulator. 
This  narrow  bandwidth  implies  that  the  optical  wave  varies  slowly  in  space  and 
time  over  one  wavelength.  This  means  we  may  assume  (E'  «  kE,  E  «  coE, 
<p'«k<p,     4>«Q)<t>),    where    (..)'  =  d(..)/dz,    and    (..)  =  d(..)/dt.        These 
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assumptions  comprise  the  slowly-varying  amplitude  and  phase  approximation 
(SVAP)  [1]. 

Neglecting  transverse  effects,  the  radiation  present  in  the  undulator  can  be 
taken  to  be  a  circularly-polarized  plane  wave  whose  vector  potential  is 

E(t) 


A  = 


[sini//,cos  y/,0] 


(2-17) 


where  w  =  kz-cot  +  <p.    Using  this  vector  potential  in  the  Coulomb  gauge  of  the 
wave  equation  yields 


(  1    ri2  ^ 

V2- 


V 


&  dt' 


JdE      1<?EV  nl 


) 


\dz     c  dt 


-2E 


(  a 


C(f)       1  o(f> 


Att-. 


^r  +  —£  \isin  v.cos  if/,0]  =  — J± 
\dz     c  dt r  c 


(2-18) 


where  the  SVAP  approximation  has  been  used  to  eliminate  all  terms  with  two 
derivatives.  The  total  beam  current  density  J2  is  the  sum  of  all  the  single 
particle  currents  J2j  =-ecpjSz'{x  -rt),  where  ?.,  is  the  trajectory  of  the  ith 
electron  and  fiu  is  given  in  equation  (2-10).  The  fast  rotating  sine  and  cosine 
functions  of  y/  can  be  eliminated  by  defining  two  orthogonal  unit  vectors 
^  =[cos  y,-sin  y/,0]  and  e2  =[sin  y/,cos  y/,0].  Projecting  (2-18)  onto  these  two 
unit  vectors  separates  the  wave  equation  into  two  first-order  scalar  differential 
equations 


and 


dE      \dE      27T-,     . 
dz      c  dt        c     ± 


Jab      *\dd>\      2n-,     . 

E  — -  +  — -    =  — J ,  •  e2 
[dz     c  dt)       c  ■ 


(2-19) 


(2-20) 
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Making  the  change  in  coordinates  z  =  z+ct  in  equations  (2-19)  and  (2-20) 
simplifies  the  differential  operator  to  \dj dz  +  c~^  d,  di)  ->  c'^a/at.  The  new 
coordinate  z  stays  with  a  point  on  the  field  envelope  traveling  at  speed  c.  The 
electrons  then  drift  back  with  a  relative  speed  c(l-/?2).  Inserting  this  change,  as 
well  as  the  expression  for  the  current  density,  yields 


6E  =     2;zeKc 

at  '        y 

and 


cosf^+v)^31  (*-/;)  (2-21) 


E^=2^Kcsin(y/+       }^(._?j)  (2_22) 

at  y 

Next  we  substitute  y/  +  ^  =  £+  <t>- 

A  volume  element  much  smaller  than  the  coherence  volume,  but  larger  than 

one  wavelength  of  light,  is  selected  to  sample  electrons.    Since  the  fields  are 

slowly  varying,  only  one  wavelength  of  phase  space  is  selected  and  sampled 

from  this  volume  element,  because  all  wavelengths  in  the  volume  must  evolve  in 

the  same  way  [1].  The  electrons  in  one  wavelength  of  phase  space  are  sampled 

to  achieve  an  accurate  representation  of  initial  phase  space  distribution,  so  the 

sum  over  all  particles  in  equations  (2-21 );  and  (2-22)  can  be  replaced  by  an 

average  over  the  sampled  electrons  weighted  by  the  electron  density  p  [1]. 

Indicating  the  ensemble  average  of  particles  at  a  site  z  by  (..),  the  optical  wave 

equations  become 


6E  _  IneKcp 
at  "         y 

and 

a<fr  2  7ieKcp 


cos 


ti+  4))  (2-23) 


E2BfS!£(^+^.  (2-24) 

ot  y       x  ' 


14 


Finally,  by  introducing  the  dimensionless  time  r  =  ct/L,  current  density 
j  =  8N(e7iKL)2p/ylmc2 ,  and  field  strength  \a\  =  A7rNeEKL/y20mc2,  equations  (2- 
23)  and  (2-24)  further  simplify  to 

\a\=  -j(cosU+  <f)),  c</>=±!sin(C+  <fi).  (2-25) 

In  low  efficiency  FELs,  the  energy  loss  per  electron  within  the  undulator  is  small, 
and  so  we  hold  the  electron's  energy  constant  at  its  initial  value  of  y0.   In 
complex  notation,  the  wave  equation  is 

a  =  -y(e-*),  (2-26) 

where  a  =  \a\e  [1].  Equations  (2-25)  and  (2-26)  show  that  electron  phases  near 
£+  <f>K  7r  drive  the  optical  amplitude  and  lead  to  gain,  while  the  optical  phase  is 
driven  by  electron  phases  near  <^+0«;r/2.  These  equations  complete  the 
feedback  mechanism  of  the  FEL.  The  pendulum  equation  describes  the 
bunching  of  electrons  which  drives  the  optical  amplitude,  and  the  optical 
amplitude  controls  the  amount  of  bunching  that  takes  place.  For  the  linear 
undulator,  the  dimensionless  current  j  is  modified  by  letting 
K  ->  K(j0(<Z)  -  J,{£))  where  J0  and  J,  are  Bessel  functions  and  £=  K2/2(l  +  K2). 

The  previous  derivation  of  the  wave  equation  disregarded  the  transverse 
evolution  of  the  fields  by  making  the  plane-wave  approximation.  To  include 
these  transverse  diffraction  effects,  the  three-dimensional  FEL  wave  equation 
uses  the  SVAP  approximation  in  a  similar  manner  to  equations  (2-26)  and 
results  in  the  parabolic  wave  equation  [8] 

(^Vi+|)a(*,r)  =  >-'),  (2-27) 
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where  vJ  =  ^  +  <^,  and  transverse  coordinates  are  normalized  to  the 
characteristic  mode  size,  so  that  x(x/LX)v2  ->x,  and  y(^/LA)V2  -»y,  and  the 
longitudinal  coordinate  is  (z-ct)/NA->z.  Equation  (2-27)  can  be  used  to 
include  the  effects  of  electron  beam  radius  and  diffraction  effects. 


D.  FEL  DIMENSIONLESS  PARAMETERS 

Dimensionless  parameters  are  extremely  useful  when  discussing  FEL 
attributes.  Several  of  them  have  already  been  defined  and  used,  such  as  the 
undulator  parameter  K,  the  electron  phase  and  phase  velocity  Q  and  v,  the 
current  density;,  and  the  optical  field  a. 

The  general  method  of  normalizing  consists  of  finding  characteristic  lengths 
and  times,  used  as  divisors  to  remove  all  dimensions.  The  characteristic  time 
scale  is  given  by  the  time  for  an  electron  to  traverse  the  undulator,  L/c.  The 
characteristic  longitudinal  length  for  macroscopic  processes  is  the  undulator 
length,  L.  The  characteristic  transverse  macroscopic  length  is  the  optical  mode 
waist  size  [1].  The  characteristic  microscopic  longitudinal  distance  is  the 
slippage  distance,  or  the  distance  that  light  passes  over  an  electron  bunch  as  it 
travels  down  the  undulator.  Looking  at  an  electron  -  photon  race  down  the 
undulator,  the  photon  beats  the  electron  by  a  distance 
As  =  Av-Ar  =  (c-/?zc)L/c  =  (l-/?z)/V/i0.  Using  (1-/?Z)  *  (l  +  K2)/2/2,  for 
relativistic  electrons  combined  with  the  resonance  condition  (2-16),  the  slippage 
distance  can  be  written  as  As  »  NA.  The  slippage  distance  is  the  characteristic 
length  over  which  an  electron  and  photon  can  exchange  energy. 

Finally,  the  Rayleigh  length,  Zq,  is  a  measure  of  the  optical  beam  diffraction 
[1]  and  is  determined  by  the  resonator  configuration  including  mirror  separation 
and  curvature.    The  optical  mode  waist  is  related  to  the  Rayleigh  length  by 
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ttW*  =xZ(j.  The  dimensionless  Rayleigh  length  is  given  by  z0  =Z0/L  [1],  and  the 
optical  mode  waist  is  normalized  by  JEkfn  so  that  vv0  =Jz^- 
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III.  QUANTUM  LIMITATION  ON  FREE  ELECTRON 
LASER  START-UP 


A.  INTRODUCTION 

A  natural  place  to  begin  examining  the  FEL  is  at  start-up.  This  chapter  looks 
into  processes  that  occur  during  the  first  few  passes  of  the  electrons  through  the 
undulator.  For  special  FEL  designs,  the  number  of  photons  in  the  resonator  may 
be  few  enough  to  cause  quantum  mechanical  effects  in  this  otherwise  classical 
laser. 

On  each  pass  through  the  FEL  oscillator,  the  periodic  magnetic  fields  of  the 
undulator  give  the  electron  a  transverse  acceleration,  causing  it  to 
spontaneously  emit  radiation  into  the  resonator  modes  [9].  On  subsequent 
passes  through  the  oscillator,  the  electrons  continue  to  spontaneously  emit 
radiation,  but  now  in  the  presence  of  existing  radiation.  This  results  in 
stimulated  emission  and  absorption  of  photons  by  the  electrons.  As  the  number 
of  photons  in  the  resonator  increases  sufficiently  to  describe  a  classical  field, 
the  radiation  frequencies  with  the  largest  gain  grow  to  dominate  the  resonator 
mode  and  the  spectrum  narrows.  If  the  single-pass  gain  of  the  FEL  is  above  the 
resonator  losses  (the  "threshold  condition"),  the  radiation  power  grows 
exponentially  to  saturation  in  strong  fields.  For  most  operating  FELs,  the  start- 
up process  occurs  in  a  few  hundred  passes  However,  a  new  class  of  compact 
FELs  with  short  undulator  periods  Z0,  a  large  number  of  periods  N ,  and  a  low 
undulator  parameter  K  =  eBA0/27rmc2 ,  may  experience  significant  delays  in 
start-up  due  to  large  fluctuations  caused  by  the  low  number  of  photons 
spontaneously  emitted  into  the  resonator  mode  each  pass.    In  some  cases,   the 


FEL  reaches  steady-state  without  sufficient  photons  to  achieve  a  classical  field, 
and  the  FEL  fails  to  operate. 


B.  COMPACT  FEL  DESCRIPTION 

Compact  FELs  are  a  new,  exciting  direction  in  FEL  research  and  offer 
attractive  advantages  over  conventional  FELs  in  several  ways.  The  compact 
FEL  has  a  significantly  reduced  undulator  period,  A0  a  1mm,  which  should 
reduce  the  size,  weight,  and  cost  of  the  FEL.  The  shorter  undulator  period  also 
allows  production  of  a  short  wavelength  radiation  with  a  less  powerful 
accelerator,  since  the  radiation  wavelength  is  inversely  proportional  to  the 
electron  beam  energy  ymc2,  X  =  AQ^  +  K2)/2y2 .  Unfortunately,  decreasing 
the  undulator  period  reduces  the  coupling  of  the  electrons  to  the  transverse 
radiation  field,  since  the  deflection  of  the  electrons  in  the  undulator  is  given  by 
Ax  *  KAQ/27iy.  To  make  matters  worse,  there  are  technological  problems  with 
obtaining  high  magnetic  field  densities  in  very  short  undulator  periods,  thereby 
reducing  the  undulator  parameter  to  a  value  K  «1.  This  problem  acts  to  further 
reduce  the  gain.  The  length  L  =  Naq  of  the  compact  undulator  must  be 
increased  to  compensate  for  the  reduced  gain. 

The  accelerator  used  by  the  compact  FEL  to  produce  short  wavelength 
radiation  may  be  much  less  energetic  than  that  required  by  a  typical  FEL  since 
the  undulator  period  is  reduced.  But,  because  the  compact  FEL  has  an 
increased  number  of  periods,  the  accelerator  must  have  higher  beam  quality 
consistent  with  the  decreased  gain  spectrum  bandwidth,  1/(2A/).  In  this 
analysis,  the  beam  quality  is  considered  achievable,  so  that  neither  emittance 
nor  energy  spread  degrades  the  FEL  gain.  In  this  compact  FEL  design,  the 
undulator  period  is  small  which  leads  to  a  small  undulator  parameter  K  «1 .   A 
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large  number  of  undulator  periods  are  used  to  insure  the  expected  classical  gain 
is  above  the  threshold  condition. 


C.  START-UP  OF  THE  COMPACT  FEL 

On  each  pass  through  the  undulator,  the  number  of  spontaneously  emitted 
photons  from  a  single  electron  into  the  resonator  mode  is  calculated  as  follows. 
The  total  power  radiated  by  a  single  electron  is  given  by  [10] 

2eVtf2 
P  =  ze  r  r±  (3-1) 

3c 
where  fix  =  Kk0c/y  is  the  electron's  transverse  acceleration.    The  energy  of  a 

photon  emitted  at  the  radiation  frequency  satisfying  the  resonance  condition, 

a>  =  2y2k0c,  is 

E*2hck0y2.  (3-2) 

The  number  of  photons  emitted  by  one  electron  in  one  pass  through  the 
undulator  is  then  obtained  by  multiplying  the  power  equation  (3-1)  by  the  time  an 
electron  spends  in  the  undulator,  Af«L/c,  and  then  dividing  by  the  photon 
energy  to  obtain 

WT  =  2NaK2.  (3-3) 

Here  a  =  e2/hc  is  the  fine  structure  constant. 

The  quantity  in  equation  (3-3)  is  the  number  of  photons  emitted  into  the  solid 
angle  y~2  [10].  As  Figure  3-1  shows,  this  solid  angle  is  typically  large  enough 
so  that  most  photons  miss  the  mirror  at  the  ends  of  the  resonator  and  are  lost. 
In  order  to  calculate  the  fraction  of  the  photons  captured  in  the  resonator,  we 
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must  multiply  WT    by  the  ratio  of  the  solid  angle  subtended  by  the  fundamental 
resonator  mode. 


Figure  3-1 .  Diagram  showing  relative  sizes  of  the  emission  angle  y"1  and  the  largest  angle  at 
which  photons  can  be  stored  within  the  resonator,  ec. 


The  optical  mode  radius  is  given  by  w2  =w02(l+z2/z2,),  where  vv0  is  the  optical 
mode  waist,  z0  =  7iw2  /  /,  is  the  Rayleigh  length,  so  that  the  solid  angle 
subtended- by  the  optical  mode  is  «  l/(2/v>2).  Multiplying  the  number  of  photons 
emitted  by  an  electron  per  pass,  WT,  by  the  ratio  of  solid  angles  gives  Wr  =  aK2 
[3]  as  the  number  of  photons  emitted  per  pass  that  are  saved  in  the  resonator 
mode. 

Each  electron  in  the  beam  emits  WY  photons  into  the  resonator  mode.  The 
volume,  Vc  *  N2A2A0,  is  used  to  count  electrons  and  photons  and  is  determined 
by  a  longitudinal  length  corresponding  to  the  gain  bandwidth  NX,  and  the 
transverse  mode  area  La=NSJ.0  [1],  The  number  of  photons  spontaneously 
emitted  in  one  pass  into  this  volume  is  Sr  =  pFaK2N2?2?.Q  where  p  is  the  electron 
beam  particle  density,  F  =  ttt2 ^/^w02  is  the  filling  factor  describing  the  transverse 
overlap  of  the  electron  and  optical  beams,  rb  is  the  electron  beam  radius,  and 
tv0  «(LA/3;r)v   is  the  radiation  mode  waist . 
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After  a  few  passes,  there  are  Ny  photons  in  the  volume  Ve.  The  resulting 
radiation  field  follows  a  Poisson  distribution  giving  a  fluctuation  in  the  photon 
number  of  N~v2  [11].  If  these  fluctuations  are  larger  than  the  gain  bandwidth  of 
the  FEL,  1/(2A/),  they  act  to  reduce  gain.  To  include  the  effect  of  these  photon 
fluctuations  we  consider  a  simple  analogy.  Since  the  dimensionless  radiation 
field  is  given  by  \a\  -  A7rNeKLE/ylmc2 ,  where  E  is  the  radiation's  electric  field,  a 
fractional  fluctuation  in  E  is  equivalent  to  a  fractional  fluctuation  in  the  undulator 
parameter  K .  Thus,  to  the  electron,  the  radiation  fluctuations  are  equivalent  to 
an  imperfect  undulator  field  of  A/J/2  coherent  periods.  Clearly,  A/r1/2  «  N 
signals  the  presence  of  a  classical  field  within  the  resonator.  The  gain  in  the 
undulator  is  proportional  to  the  dimensionless  current  density 
j  =  8N3(exKA0)2pF/ylmc2  ,  although  j  must  be  corrected  by  the  ratio 
[n\/2 /n)  for  the  calculated  gain  to  account  only  for  the  coherent  number  of 
periods  in  the  undulator.  However,  there  are  [N/Nry2J  coherent  undulators 
acting  in  series.  Thus,  the  actual  gain  per  pass  after  one  pass,  including 
quantum  fluctuations,  is  estimated  as 


G  =  G0T^,  (3-4) 


where  G0  «  0.135)  is  the  classical  small-signal  gain.  If  the  number  of  photons 
spontaneously  emitted  on  the  first  few  passes  is  low,  the  photon  fluctuations  may 
be  large  and  lead  to  a  significant  decrease  in  gain. 

After  many  passes,  the  number  of  photons  in  the  resonator  will  change  due 
to  continued  spontaneous  emission,  the  resonator  Q  losses,  and  possible  gain 
from  stimulated  emission.  The  differential  equation  describing  the  evolution  of 
the  photon  number,  Nr,  over  many  passes,  n,  is  [11] 
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^  =  S,  +  Wr(G-Q-1),  (3-5) 

where  Q  represents  the  loss  within  the  resonator.  The  FEL  may  remain  below 
threshold,  G  <  Q"1,  and  yield  a  steady  state  when  the  number  of  photons 
reaches 
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(3-6) 


Equation  (3-6)  shows  the  possibility  of  a  complex  number  of  photons,  which  can 
be  interpreted  as  the  absence  of  a  steady-state  solution  and  no  start-up  of  the 
FEL.  From  this  interpretation,  it  is  inferred  that  the  start-up  condition  is 
4SyG0Q2/N2  >  1.  Substituting  Sr  -  pFaK2N2??A0,  the  dimensionless  current 
density  /,  the  small  signal  gain  G0~0.135/,  and  introducing  the  Compton 
wavelength  Xc  =  h/mc   leads  to  a  general  start-up  condition 

0.135(yQ)2^  >  {2tt)2N3Xc.  (3-7) 

In  the  classical  limit,  Xc  -^0,  all  FELs  satisfy  this  requirement  and  start.  A 
typical  FEL  (such  as  the  original  Stanford  SCA/FEL  designed  with  jQ*(4x)2, 
Lorentz  factor  y  ^80,  A  =  3{im,  and  N  =  160)  satisfies  the  requirement  of 
equation  (3-7)  by  two  orders  of  magnitude. 

On  the  other  hand,  a  compact  FEL  with  parameters  of  y&24,  A  =  86^m, 
jQ -70,  and  A/  =  103  results  in  no  start-up.  Two  figures,  3-2  and  3-3,  appear 
below.  Figure  3-2  represents  the  results  of  a  conventional  ("classical") 
calculation  which  does  not  include  quantum  fluctuations.  This  calculation  and 
that  figure  erroneously  predict  successful  start-up.  Figure  2-3,  which  does 
include  quantum  fluctuations,  correctly  predicts  that  the  FEL  does  not  start  up. 
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The  numerical  simulation  including  the  classical  spontaneous  emission  line 
spectrum  is  used  to  determine  the  photon  evolution  per  pass  without  the 
quantum  fluctuation  term  included  in  the  spectrum.  The  bottom  panel  shows  the 
classical  spontaneous  emission  spectrum  while  the  middle  panel  shows  the 
photon  number  evolution  over  n=300  passes  with  the  final  spectrum  at  the  top. 
Over  a  few  passes,  the  photon  number  spectrum  shifts  to  the  maximum  gain 
frequency  and  narrows  to  a  classical  FEL  line  spectrum.  The  gain  spectrum  Gy 
is  displayed  with  the  dotted  line  showing  the  loss  per  pass.  The  total  number  of 
photons  spontaneously  emitted  per  pass  is  given  by  Sv  =278.  The  total  number 
of  photons  in  the  final  laser  spectrum  is  N  =  2.9  x  1015.  Normally,  high-power 
saturation  would  limit  the  final  number  of  photons  to  a  value  less  than  N^  before 
n  =  300  passes,  but  high-power  saturation  is  not  included  in  this  simulation. 
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Figure  3-2.  Photon  number  evolution  using  small-signal  gain. 
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Figure  3-3  displays  the  results  of  including  quantum  fluctuations  in  the  gain 
term.  After  n-2000  passes,  the  upper  panel  shows  the  photon  number  spectrum 
begins  to  increase  in  amplitude  and  evolve  toward  the  maximum  gain  frequency. 
Before  reaching  the  threshold,  a  steady  state  is  reached  with  Nrf  =  2.7  x104 
photons,  where  the  increase  in  photons  per  pass  due  to  spontaneous  emission 
Sv  =278,  and  stimulated  emission  GNr,  is  equal  to  the  losses  in  the  system, 
QNy.  The  actual  gain  of  the  laser,  including  the  quantum  fluctuation  term,  is 
below  the  losses  of  the  laser,  as  shown  in  the  center  panel.  Although  this 
compact  FEL  would  be  expected  to  work  according  to  classical  FEL  theory, 
Figure  3-3  shows  a  steady  state  photon  number  is  reached  in  the  resonator 
below  that  required  for  a  classical  field.  Thus,  the  gain  never  exceeds  the 
losses  and  the  laser  does  not  start. 
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Figure  3-3.  Photon  number  evolution  using  the  modified  equation  (3-4)  as  the  gain. 
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D.  ANALYTICAL  RESULTS 

Assuming  a  form  for  the  gain  of  G  *  GQNr/N2,  the  photon  rate  equation  (3- 
5)  can  be  solved  analytically  for  the  photon  number  evolution  per  pass.  Again, 
this  solution  is  valid  only  up  to  Nr  =  N2  at  which  time  the  model  used  to  develop 
the  gain  formula  is  no  longer  valid.  The  solution  separates  into  two  cases:  one 
when  the  start-up  condition  (3-7)  is  satisfied,  and  the  other  when  it  is  not.  If  the 
start-up  condition  is  not  satisfied,  the  photon  number  after  the  nm  pass  is  given 
by 
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(3-8) 


where  q  =  (1  -  agqsyq7 /n2)':  and  r  =  |(1  +  q)/(q - 1)|.  This  solution  displays  a 
linear  growth  dependence  on  n  for  small  n.  However,  there  are  never  enough 
photons  in  the  resonator  for  the  gain  to  overcome  FEL  loss  defined  by  Q,  so  a 
steady-state  is  reached  with  an  insufficient  number  of  photons  in  the  resonator 
to  form  a  classical  field.  If  the  start-up  condition  (3-7)  is  satisfied,  the  photon 
number  after  thenm  pass  is 
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where  p  =  (4G0SVQ2//V2  - 1)  .  This  solution  for  the  photon  number  evolution 
grows  slowly  per  pass  initially  due  to  contributions  primarily  from  spontaneous 
emission.  As  the  number  of  photons  in  the  resonator  grows,  stimulated  emission 
begins  to  add  photons  and  their  number  increases  dramatically.   This  continues 
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until  the  fluctuations  in  photon  number  are  much  less  than  the  gain  spectrum 
bandwidth  resulting  in  classical  gain. 

The  threshold  for  which  stimulated  emission  produces  as  many  photons  as 
does  spontaneous  emission  is  obtained  by  setting  Sr  =NrG.  Dattoli  [12]  gives 
an  expression  for  this  threshold  number  of  photons  as  (Nr)  =  y^/[32^NXG) 
using  a  quantum  mechanical  formalism.  This  expression  can  be  reproduced 
with  all  the  correct  functional  proportionalities  much  more  simply  by  setting  the 
gain  equal  to  the  small  signal  gain  G  =  0.135).  However,  when  quantum 
fluctuations  are  included,  and  the  gain  is  determined  by  (3-4),  the  condition  for 
the  stimulated  emission  to  equal  the  spontaneous  emission  is 


(N,)    . 
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This  new  condition  is  different  from  Dattoli's  and  takes  into  account  the  reduced 
gain  initially  present  due  to  fluctuations  in  the  radiation  field.  It  can  be  seen  that 
if  the  laser  start-up  condition  is  not  fulfilled,  the  steady-state  photon  number 
within  the  resonator,  A/*,  will  be  less  than  the  threshold  number  of  photons 
required  for  equal  spontaneous  and  stimulated  emission  rates. 


E.  CONCLUSIONS  AND  FURTHER  RESEARCH  SUGGESTIONS 

The  type  of  FEL  that  may  exhibit  start-up  problems  due  to  quantum 
fluctuations  is  one  that  attempts  to  produce  short  optical  wavelengths  with  a  low 
energy  electron  beam.  The  low  energy  electron  beam  necessitates  a  small 
period  undulator,  which  in  turn  reduces  the  undulator  parameter.  Both  of  these 
properties  act  to  reduce  the  expected  gain  of  the  FEL.     In  order  to  achieve 
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sufficient  gain,  the  number  of  periods  in  the  undulator  must  be  increased. 
However,  this  design  logic  can  lead  to  an  FEL  that  may  not  operate  due  to 
quantum  fluctuations  during  start-up. 
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IV.  ANALYTICAL  DESCRIPTION  OF  LONGITUDINAL 

MODES 


A.  INTRODUCTION 

The  optical  power  in  the  FEL  oscillator  starts  from  spontaneous  emission 
and  grows  to  strong-field  saturation  over  many  passes.  During  the  crucial 
startup  process,  coherence  develops  in  the  weak  optical  field  and  small 
microscopic  bunching  develops  in  the  electron  beam.  In  this  chapter  the  FELs 
considered  are  typical  in  the  sense  that  the  quantum  effects  previously  explored 
do  not  occur.  The  FELs  considered  here  have  weak  classical  fields,  low  gain, 
and  electron  pulses  much  longer  than  the  slippage  distance  NA. 

Assuming  low  gain,  a  long  pulse,  and  weak  optical  fields,  Elleaume  [12] 
showed  that  FEL  optical  pulse  evolution  could  be  described  by  a  Schroedinger- 
type  equation  with  a  non-Hermitian  harmonic  oscillator  Hamiltonian  .  Using  the 
same  assumptions,  Dattoli  [13.14,15,16]  employed  methods  of  quantum  optics  and 
an  integro-differential  equation  to  obtain  his  solution  to  the  same  problem.  In 
each  case,  the  mathematical  result  involves  second-order  derivatives.  In  this 
chapter,  we  expand  the  optical  and  electron  pulse  shapes  in  a  Taylor  series  to 
derive  a  first-order  partial-differential  equation  describing  the  optical  field 
evolution  over  many  passes  through  the  undulator.  This  equation  is  used  to 
investigate  coherence  development  and  the  effects  of  the  electron  pulse  shape 
on  the  optical  field  amplitude  and  phase. 
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B.  MATHEMATICAL  DEVELOPMENT 

In  order  to  generalize  the  optical  field  in  the  longitudinal  direction,  to  allow 
for  multiple  optical  modes,  we  must  follow  the  optical  field  envelope  evolution  all 
along  z.  The  extension  to  multi-mode  theory  by  transforming  a  ->  a{z)  is 
equivalent  to  allowing  the  optical  field  a  continuum  of  wavenumbers  a  -»  a(/c). 
Fully  extending  the  theory  uses  a-+a(z),  v->  v(z),  and  £->£(z).[i]  As 
discussed  in  Chapter  II,  all  longitudinal  distances  are  normalized  to  the  slippage 
distance,  the  distance  over  which  the  light  and  electrons  can  exchange 
information,  so  that  z/NA  ->  z.  The  electrons  slip  back  one  slippage  distance, 
or  Az  =  1,  with  respect  to  the  light  during  the  interaction  time  r=0  ->1  leading  to 
NX  =  Az  =  r. 

The  FEL  interaction  can  be  described  by  the  optical  wave  equation  [1], 

0 

a(z,  T)  =  jo(z  +  r)(exp[-/^z  +  r,  r)])  (4-1 ) 

coupled  to  the  electron  pendulum  equation  [1], 

oo  o 

£(z,  t)  =  v(z,  r)  =  |a(z  -  r,  t)\cos[£z,  t)  +  <f>{z  -  t,  r)],  (4-2) 

o 

where      (..)  =  d(..)/ dr,      £=  (k  +  k0)z-  cot       is      the      electron      phase, 

o 

v-  L[(k  +  k0)pz  -  k]  =  C,  is  the  electron  phase  velocity,  |a(z,  r)|  is  the  optical 
electric  field  amplitude,  4{z,z)  is  the  optical  phase,  and  <..>  represents  an 
average  over  all  of  the  electrons  in  the  beam.  The  electrons'  pulse  shape  is 
described  by  o(z)  with  peak  value  of  unity.  The  longitudinal  field  sites  z  refer  to 
positions  on  the  optical  field  envelope.  In  (4-1 ),  the  light  remains  fixed  in  z  while 
the  electrons  slip  back  to  a  field  site  z  +  r,  while  in  (4-2),  the  electrons  remain 
fixed  and  the  light  slips  ahead  by  z  -  r  as  it  moves  through  the  undulator.   This 
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slippage  allows  the  light  and  electrons  to  exchange  information,  and  allows  the 
electrons  to  transfer  of  information  from  one  optical  site  to  another. [1] 

When  the  FEL  is  operating  in  the  low-gain  regime,  we  expand  equation  (4-1) 
to  lowest  order  in  the  current  density,  j,  and  assume  the  field  changes  are  small 

O  O 

over  a  single  pass,  a|«1a|  and  (f>«<j>.  Since  the  optical  fields  are  weak,  the 
electron  phase  may  be  expanded  to  lowest  order  in  the  field  magnitude  |a|, 
giving  £=£0  +  v0z+...  .  Here  £0  is  the  initial  electron  phase,  and  v0  is  the  initial 
electron  phase  velocity.  The  long  pulse  shapes  may  be  expanded  in  a  Taylor 
series  in  the  normalized  slippage  distance,  r:  o(z+  r)=  o(z)+  ro/(z)+..., 
|a(z- r)|  =  |a(z)|- zja(z)|'+...  and  <p(z-  r)  =  </>(z)  -  r<p'{z)+...  where 
(..)'=  d(..)/dz.  These  assumptions  allow  the  expansion  and  integration  of 
equation  (4-2), 


\d\ 

&z  +  r,  t)  =  £0  +  vQt  +  J-^-[-cos(4"0  +  v0t+  <p)  +  cos(4o  +  <f>)  -  v0TSin{£0  +  </>)] 

vo 

[2sin{£0  +  v0r+  <j>)  -  2sin{£0  +  <p)-2  vQrcos{£0  +  <f>)  v0V(C0  +  j)] 


\a\<p' 


+  i-ly-  [2  cos(  £0+v0r +</>)- 2  cos{  £0  +  <f>)-2v0  rsin{  &  +  <p)  +  v02  r2  cos{  £,  +  ^)] 
vo 

+ ,  (4-3) 

where  |a|  =  |a(z)|,  \a\   =  |a(z)|  ,  <f>=  <p(z),  and  ft  -  <f>'{z). 

For  a  monoenergetic  beam  (where  all  electrons  have  the  same  initial  phase 
velocity,    v0)   the   electron   dynamics   expressed    in   equation   (4-3)    can    be 


incorporated  into  the  optical  field  equation  (4-1)  using  the  phase  average 
Jo 


(..)  =  J     (..)d£0  /2tt.  The  complex  optical  field  equation  then  becomes 
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where  a  =  |a(z)|exp(/'^(z)).  This  equation  includes  the  electron  dynamics  of 
equation  (4-2)  using  the  assumption  of  low  gain,  weak  fields,  and  long  pulses  as 
expressed  in  equation  (4-3). 

Integrating   over  a   single   pass,    z  =  0->1,   yields  the  following   partial- 
differential  equation  for  the  net  change,  Aa,  in  the  complex  optical  field, 


Aa  =  jaog(v0)-±-ja-?-h(v0)--!-ja—g,(v0)         (4-5) 
2      oz  2      cz 


with  constant  coefficients  for  a  given  initial  electron  phase  velocity 


9  = 


2  -  2e-'v°  -  iV0(l  +  e"Vb )        6  -  6e_'v'0  -  6/V0e-/V°  +  vg(l  +  2e-,v° ) 


2£ 


,h  = 


,_  * 


2vn 


^'  = 


CfVn 


(4-6) 


With  no  pulses,  a  =  1  and  a' =0  so  there  are  no  longitudinal  modes.  The 
second  and  third  terms  on  the  right  side  of  equation  (4-5)  are  zero.  In  that  case, 
equation  (4-5)  describes  the  usual  single-mode,  low-gain  evolution  of  the 
complex  optical  field  [1].  The  real  part  of  g(v0)  is  the  familiar  antisymmetric  gain 
spectrum  in  v0,  and  the  imaginary  part  describes  the  symmetric  shift  in  the 
optical  phase  [1].  The  phase  velocity  v0  determines  the  resonance  between  the 
optical  wavelength,  undulator  wavelength,  and  electron  z  velocity;  maximum  gain 
occurs  for  v0  =2.60616.  The  functions  h(v0)  and  g'{  v0)  in  (4-5)  determine  the 
effects  of  slippage  and  the  electron  pulse  shape  on  the  optical  field  evolution. 
Elleaume  [13]  and  Dattoli  [14,15,16,17]  obtained  similar  functions  of  v0  in  their 
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treatment  of  the  long  pulse  problem,  but  they  also  keep  some  second-order 
derivatives  in  the  expansion  of  the  field. 

Over  many  passes  n  through  the  undulator,  the  effects  of  resonator  losses 
and  desynchronism  must  be  included  [1].  Desynchronism  measures  the 
mismatch  at  the  beginning  of  the  undulator  between  the  arrival  times  of  the 
succession  of  electron  pulses  from  the  accelerator,  and  the  rebounding  optical 
pulse  in  the  resonator  cavity.  The  dimensionless  desynchronism  is 
d  =  -2 AS/ NPl,  where  AS  is  the  difference  between  the  actual  resonator  length 
and  the  synchronous  cavity  length  [1].  The  resonator  loss  over  many  passes  is 
described  by  the  resonator  quality  factor  Q.  The  resulting  first-order  partial 
differential  equation  gives  the  development  of  longitudinal  modes  in  a  long  pulse 
over  many  passes  n  in  weak  fields  of  a  low  gain  FEL, 

da      .       ,     .     /  .  do ..     .     /  .   da    , ,     ,      a       .3a     ..-.. 
^  =  ya^0;--,a-/,K;--ya-gK;---d-.    (4-7) 

The  complicated  optical  pulse  evolution  determined  by  equation  (4-7)  depends 
only  on  the  dimensionless  current  /  the  resonator  Q,  the  desynchronism  d,  the 
electron  pulse  shape  c(z),  the  initial  optical  pulse  amplitude  \a0(z)\  and  phase 
<fi0(z),  and  the  initial  phase  velocity  v0. 

Since  the  shape  of  the  phase,  ^(z),  is  free  to  evolve,  the  phase  velocity  v0  is 
somewhat  arbitrary.  For  instance,  if  the  phase  has  a  linear  slope,  <f{z)  =  Akz,  it 
is  equivalent  to  a  change  in  wave  number,  k  ->k  +  Ak,  and  alters  the  effective 
value  of  the  phase  velocity.  Therefore,  the  phase  velocity  v0  in  equation  (4-7) 
can  be  chosen  for  convenience  with  no  loss  of  generality.  A  unique  and 
important  value  of  v0  is  the  phase  velocity  giving  peak  FEL  gain,  gr'(v0)  =  0. 
The  optical  pulse  will  naturally  evolve  towards  v*  =2.60616  because  of  mode 
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competition.  In  fact,  v0  cannot  be  chosen  far  from  vj,  because  the  slope  of  the 
phase  required  to  evolve  to  peak  gain  can  exceed  the  restrictions  of  the  slowly- 
varying  amplitude  and  phase  approximation.  This  phenomenon  can  be  seen  by 
beginning  with  the  phase  velocity  v0  =  L[(k  +  k0)pz  -  k\  and  deriving  the  change 
in  phase  velocity  for  a  change  in  wave  number  Av=  -A/c  =  -$ .  It  can  be  seen 
that  the  slope  of  the  phase  provides  the  change  in  wave  number.  Next,  using 
the  slowly-varying  phase  approximation,  the  change  in  the  slope  of  the  phase 
must  be  small  as  the  optical  field  makes  one  pass  in  the  undulator,  so  AvAr<  1. 
During  one  pass,  the  optical  phase  interacts  over  a  dimensionless  distance  of 
Az  which  can  be  linked  to  the  size  of  the  electron  pulse  oz.  This  in  turn  gives  a 
new  requirement  A  vAoz  <  1.  Since  one  of  the  basic  assumptions  of  this  analysis 
that  oz  »  1,  the  restriction  Av«^  is  implied.  When  v0  =  v0  and  gr'=0,  the 
equations  for  the  amplitude  and  phase  in  (4-7)  become  uncoupled,  and  are 
given  by 


^  +  (d-±jg;a)^  =  (jgra  +  ±fto>-y2Q)\a\,  (4-8) 

an  *  dz 

&  +  (d-ljg!cr)&  =  jgla-jijhro'.  (4-9) 

on  dz 

The  equations  describing  the  evolution  of  the  optical  amplitude  and  phase  can 
now  be  solved  separately. 


C.  ELECTRON  PULSE  EFFECTS 

In  this  section,  three  different  shapes  for  the  electron  pulse  are  considered 
and  then  used  to  solve  equations  (4-8)  and  (4-9)  analytically.  The  method  of 
characteristics  [18]  will  be  used  to  solve  equations  (4-8)  and  (4-9). 
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The  left-hand-side  of  equations  (4-8)  and  (4-9)  can  be  seen  to  be  the 
directional  derivative  of  |a(z,/?)|  and  <f{z,ri)  respectively,  in  the  direction  of  the 
vector  [l,(cf  -  %jg'cr)\  The  equations 


and 


ds 

jg;° 

2 

ds 

(4-10) 


(4-1 1 ) 


determine  a  family  of  curves  z  =  z(s)  and  n  =  n(s)  whose  tangent  vector, 
[z'(s),n'(s)],  coincides  with  the  vector  [l,(cf--j/gf'a)].  Therefore,  the  derivative 
of  a(z,n)  along  these  curves  becomes 
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(4-13) 


after  using  the  chain  rule  and  equations  (4-10)  and  (4-11).  The  family  of  curves 
z  =  z(s),  n  =  n(s),  and  a{s)  -  |a(s)|e,<!5<si,  determined  from  the  solutions  of  the 
system  of  ordinary  differential  equations  (4-10)  through  (4-12)  or  (4-10), (4-11), 
and  (4-13)  are  called  the  "characteristic  curves"  of  the  partial  differential 
equations  (4-8)  or  (4-9)  [18].  The  approach  of  solving  equations  (4-8)  or  (4-9)  by 
using  the  solutions  of  equations  (4-10)  through  (4-12)  or  (4-10), (4-11),  and  (4- 
13)  is  called  the  "method  of  characteristics".    The  existence  and  uniqueness 
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theorem  assumes  certain  smoothness  conditions  on  the  functions  of  z  and  n 
which  guarantee  that  exactly  one  solution  curve  (z(s),A7(s),a(s))  of  equations  (4- 
10)  through  (4-13)  passes  through  the  initial  condition  point  (z0,nQ,a0).  Using 
the  initial  values  to  create  the  parametric  initial  solution  curve  gives  z  =  z0,  n  =  0, 
and  |a|  =  |a0(z0)|  where  |a0(z0)|  is  the  initial  optical  field  shape  and  ^=^0(z0) 
where  <f>Q{z0)  is  the  initial  optical  phase.  Equations  (4-10)  and  (4-11)  can  be 
solved  for  z  and  n  as  functions  of  s  and  z0.  These  solutions  are  then  inverted 
and  equations  (4-12)  and  (4-13)  are  solved  and  substitutions  made  to  obtain 
ja(z,/?)|  and  <fi(z,n).  The  three  specific  cases  follow. 
1.  Continuous  Electron  Beam 

For  the  special  case  of  a  continuous  electron  beam,  o(z)  =  1  for  all  z,  at 
the  phase  velocity  for  the  peak  gain,  v0,  the  optical  pulse  can  still  evolve  over 
many  passes.  The  first-order,  partial-differential  equations  (4-8)  and  (4-9)  are 
solved  using  the  method  of  characteristics. 

Solving  equations  (4-10)  through  (4-13)  yields 

z(s,z0)  =  z0  +  (d-±jg;)s,  (4-14) 

n(s,z0)  =  s,  (4-15) 


|a(s,z0)|  =  |a(z0)|exp 


[j9r~\S 


(4-16) 


and                                     4(slzQ)  =  jgls+fo{z0)  (4-17) 
which  can  be  inverted  to  solve  for  s  and  z0  in  terms  of  n  and  z.    These  results 

are  then  substituted  into  the  expression  for  \a\  and  <f>  to  yield  solutions  for  the 
optical  amplitude  and  phase, 

\a{z,n)\  =  \aQ(z0)\exp[(jgr  -1/2Q)n] ,  (4-18) 

and 

<fr{z,n)  =  jg,n  +  <t>0(z0)  ,  (4-19) 
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where  |a0(z0)|  and  <f>0{zQ)  describe  the  initial  optical  pulse,  gr  =Re(g)  is  the  real 
part  of  g(vQ),  g,=\m(g)  is  the  imaginary  part  of  g(v0),  and  the  position 
zo  =  z-(d-2J9,')n  describes  the  optical  pulse  centroid  shifting  over  n  passes 
through  the  resonator  due  to  desynchronism.  If  d  <  ijg,  ,  the  optical  pulse  shifts 
back  over  many  passes  toward  the  trailing  edge  of  the  electron  pulse.  Further,  if 
d>ijg,  ,  the  optical  pulse  shifts  forward  toward  the  front  edge  of  the  electron 
pulse.  The  shape  of  the  optical  pulse  amplitude  and  phase  remains  the  same  as 
the  initial  pulse  throughout  the  evolution.  Figure  4-1  shows  a  comparison  of  the 
analytical  solution  (4-18)  and  (4-19)  with  the  results  of  a  numerical  simulation  for 
the  same  conditions.  The  electron  pulse  is  continuous  in  length  with  all 
electrons  at  phase  velocity  v0  for  maximum  gain. 


****      pel  Pulse  Evolution   **** 
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Figure  4-1.  Optical  pulse  evolution  comparing  a  numerical  simulation  with  the  analytical  solution 

for  a  square  electron  pulse. 
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The  desynchronism  is  d  =  0.018;  the  initially  weak  optical  field  amplitude  is 
parabolic  with  a  maximum  of  a0  =0.001  and  a  full  width  of  ua  =  10;  the  initial 
phase  is  <f>Q(z)  =  0  for  all  z;  and  the  resonator  loss  is  determined  by  Q  =  100. 
The  lower-left  panel  displays  the  continuous  electron  pulse  shape  with  a=1  for 
all  z.  The  middle-left  panel  shows  the  evolution  of  the  optical  pulse  over  n  =  60 
passes  with  the  final  optical  field  amplitude  and  shape  in  the  top-left  panel.  The 
analytical  solution  gives  a  slightly  larger  amplitude  than  the  numerical 
simulation,  but  the  centroids  and  shapes  of  the  two  solutions  are  in  close 
agreement.  The  middle-right  panel  shows  the  evolution  of  the  optical  phase 
over  the  n  =  60  passes  with  the  analytical  and  numerical  solutions  after  the  final 
pass  displayed  at  the  top-right.  The  similarity  of  the  two  solutions  is  quite 
apparent  within  the  portion  of  the  window  occupied  by  the  optical  pulse  with 
significant  amplitude.  In  regions  where  the  optical  amplitude  is  small,  the  optical 
phase  becomes  meaningless. 
2.  Parabolic  Electron  Beam 

The  parabolic  electron  pulse  has  the  form  o(z)  =  1-2z2  /  az2  for 
\z\<<jz/&,  and  o(z)=0  otherwise.  The  parabolic  electron  pulse  shape  is 
often  used  in  FEL  simulations.  In  addition,  the  parabolic  equation  expresses  the 
same  z  dependence  as  the  first  two  terms  in  a  Taylor  expansion  of  any  pulse 
shape.  Using  the  parabolic  electron  pulse  shape,  equations  (4-8)  and  (4-9)  can 
be  solved  analytically  for  arbitrary  initial  optical  pulse  amplitude  |a0(z)|,  and 
phase,  fa(z)  using  the  method  of  characteristics. 

The  solution  to  equations  (4-8)  and  (4-9)  yields  results  for  two  cases 
depending  on  whether  the  optical  pulse  shifts  forward  due  to  large 
desynchronism,  d >jjgl',  or  backward  due  to  large  gain,  d  <  ijg,' . 
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In  the  first  case,  where  desynchronism  dominates  gain,  d >±jgi'  and  the 
optical  pulse  tends  to  shift  forward,  equations  (4-10)  through  (4-13)  yield 
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where  b  =d --jjg/  and 
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Substituting  equations  (4-20)  and  (4-21)  into  equations  (4-22)  and  (4-23)  gives 
the  optical  amplitude 
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\a{z,n)\  =  \a0{z0)\exp 
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and  the  optical  phase 
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where  z0  =  \gz  b  /  jg/j     tan(q),  and  q  is  redefined  as 
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In  the  second  case,  where  gain  dominates  desynchronism,  d  <jjg/  and 
the  optical  pulse  tends  to  shift  backward,  equations  (4-10)  through  (4-13)  yield 
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\a(s,z0)\  =  \a0(zQ)\exp 
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where  w  =  2s((-£>)(yg,'))     /  crz. 

Substituting  equations  (4-26)  and  (4-27)  into  equations  (4-28)  and  (4-29) 
and  redefining  w  =  2n((-b)(jgj'))     /  <jz  yields  the  optical  amplitude 
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and  the  optical  phase 
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The  upper  sign  is  used  when  \z\<-bGz/{-bjgj')  ,  and  the  lower  sign  is  used 
when  z  is  not  between  these  limits  but  is  within  the  limits  of  the  electron  pulse. 
Outside  the  limits  of  the  electron  pulse,  the  phase  remains  constant  and  the 
optical  amplitude  has  the  form 

\a(z,n)\  =  {-[-2{z-dn)2/oj)exp{-n/2Q)  (4-32) 

where  oa  is  the  width  of  the  optical  pulse. 

The  expressions  (4-24),  (4-25),  and  (4-30)  through  (4-32),  provide  an 
analytical  solution  for  an  initially  parabolic  optical  and  electron  pulse.  Figure  4-2 
shows  the  comparison  of  the  analytical  solution  with  the  results  of  a  numerical 
solution  for  the  same  conditions.  The  desynchronism  is  d  =  0.015,  the  initial 
field  amplitude  a0  =  0.001  for  all  z,  the  phase  ^0(z)  =  0  for  all  z,  the  electron 
pulse  length  is  oz  =  10  at  phase  velocity  vQ  for  maximum  gain,  and  Q  =  100. 
The  initially  long  optical  pulse  takes  on  the  apparently  Gaussian  shape  shown  at 
the  upper  left  after  n  =  60  passes. 
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****      pel  Pulse  Evolution   **** 
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Figure  4-2.  Optical  pulse  evolution  comparing  a  numerical  simulation  with  the  analytical 

solutions  for  a  parabolic  electron  pulse. 


The  analytical  solution  gives  a  slightly  greater  amplitude  than  the  numerical 
simulation  after  n  =  60  passes,  but  the  shape  and  centroids  of  the  pulses  are 
nearly  the  same.  The  difference  in  amplitude  is  probably  due  to  an 
accumulation  of  small  numerical  errors  in  the  calculation  of  gain  each  pass.  A 
difference  in  gains  of  only  0.3%  over  n  =  60  passes  would  account  for  the 
discrepancy. 

When  d  =  0,  the  second  case  is  appropriate  because  gf/'(vj)>  0  and 
equation  (4-30)  shows  that  the  field  |a(z,n)j  -*  0  as  n  ->  <x>.  This  situation  is  the 
well-known  result  that  a  pulsed  FEL  cannot  operate  at  zero  desynchronism  [1]. 
Physically,  the  centroid  of  the  optical  pulse  is  shifting  back  over  the  electron 
pulse  by  an  amount  \\g{  for  each  pass  through  the  undulator.  Eventually  there 
is  no  overlap  between  the  pulses  and  the  optical  field  decays  to  zero.  Figure  4-3 
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shows  a  gain  versus  desynchronism  curve  for  the  parabolic  electron  pulse  and 
initially  constant  optical  pulse.  The  desynchronism  for  maximum  gain  is  located 
at  d  =jjgj',  which  is  the  desynchronism  necessary  for  the  optical  pulse  centroid 
and  electron  pulse  centroid  alignment  at  the  beginning  of  the  undulator  each 
pass. 
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Figure  4-3.  Gain  vs.  desynchronism  for  a  parabolic  electron  pulse. 


3.  Triangular  Electron  Pulse 

Another  electron  pulse  shape  that  can  be  solved  analytically  using  the 
method  of  characteristics  with  equations  (4-8)  and  (4-9)  is  the  triangular  pulse 
shape  of  the  form  o(z)  =1-|z|/cr2  with  width  oz  for  \z\  <  oz  and  o(z)=0 
otherwise.  This  case  compares  the  differences  in  optical  pulse  shape  and 
phase  caused  by  the  z2  dependence  of  the  parabolic  electron  beam  with  the  z 
dependence  of  the  triangular  electron  beam.  Due  to  the  discontinuity  in  the 
triangular  pulse  at  its  peak,  z  =0,  the  solution  has  slight  differences  in  form  for 
z  <0  and  z  >0.  An  error  near  z  =  0  is  to  be  expected,  since  the  triangular  pulse 
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shape  violates  the  assumption  that  the  electron  pulse  shape  varies  slowly; 
<?{z)  «  1.  As  in  the  parabolic  case,  the  optical  pulse  slipping  backward  due  to 
large  desynchronism,  d>^jgi\  or  forward  due  to  predominant  gain,  d >^jg' , 
appears.  However,  the  solutions  to  equations  (4-8)  and  (4-9)  for  the  triangular 
electron  pulse  shape  do  not  lend  themselves  to  being  displayed  in  a  form  that 
makes  these  two  cases  obvious.  Therefore,  we  shall  solve  for  cases  of 
\z\  <  2faz  /  jg, '  and  \z\  >  2faz  /  jg, ' . 

Solving  equations  (4-10)  through  (4-13)  for  the  first  case  of  \z\  <  2faz  /  jg,' 
yields 
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where  f  =±jg'-d.    Substituting  equations  (4-33)  and  (4-34)  into  equations  (4- 
35)  and  (4-36)  gives  the  optical  amplitude 
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The  second  set  of  analytical  solutions  are  for  \z\  >  2faz  /  /g,'.  Solving 
equations  (4-1 0)  through  (4-1 3)  for  the  case  of  \z\  >  2faz  /  jg, '  gives 
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Substituting  equations  (4-39)  and  (4-40)  into  (4-41)  and  (4^42)  gives  the  optical 
amplitude 
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where  |a0(z0)|  is  the  arbitrary  initial  optical  pulse  amplitude,  and  the  optical 
phase 
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where  <p0(zQ)  is  the  initial  optical  phase  and 
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For  both  the  amplitude  and  the  phase  solutions,  the  upper  sign  is  used  when 
z  >0  and  the  lower  when  z  <0. 

The  expressions  (4-37),  (4-38),  (4-43)  and  (4-44)  provide  an  analytical 
solution  to  the  optical  pulse  amplitude,  shape,  and  phase  in  the  presence  of  a 
triangular  electron  pulse.  Figure  4-4  shows  a  comparison  of  a  numerical 
simulation  and  these  analytical  solutions  for  the  same  initial  conditions.  The 
desynchronism  is  d  =  0,  the  initial  field  amplitude  is  a0  =  0.001  for  all  z,  the 
phase  is  <pQ(z)  -  0  for  all  z,  the  electron  pulse  length  is  oz  -  10  at  phase  velocity 
v*  for  maximum  gain,  and  Q  =  100.  The  lower-left  panel  shows  the  triangular 
electron  pulse  at  the  beginning  of  the  undulator  when  r  =  0,  and  at  the  end  when 
r  =  1  ,  showing  the  slippage  distance.    The  middle-left  panel  shows  the  optical 
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pulse  evolution  over  n  =  40  passes  with  the  final  pass  displayed  at  the  top-left. 
The  analytical  results  are  slightly  larger  in  amplitude  and  more  sharply  peaked 
than  the  numerical  simulation,  but  the  centroid  and  shape  correspond  well.  The 
middle-right  panel  shows  the  optical  phase  evolution  over  n  =  A0  passes  with 
the  final  phase  displayed  at  the  top-right.  As  can  be  seen,  there  is  excellent 
agreement  between  the  analytical  and  numerical  solutions. 


****   fel  Pulse  Evolution  **** 
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Figure  4-4.  Optical  pulse  evolution  comparing  a  numerical  simulation  with  the  analytical 

solutions  for  a  triangular  electron  pulse. 


D.  COHERENCE  DEVELOPMENT 

The  evolution  of  the  optical  amplitude  and  phase  over  many  passes  can  be 
used  to  analytically  examine  coherence  development.     Taking  the  parabolic 
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electron  beam  with  d<  i/g,'  as  an  example,  the  optical  amplitude  solution  (4-30) 
on  the  first  pass  is  given  by  |a(z,n)|  =  |a0(z0)|  where  |a0(z0)|  represents  an 
arbitrary  optical  pulse  shape  as  a  function  of  z0{z).  After  many  passes, n-»  oo, 
the  expression  for  z0  becomes 
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To  obtain  (4-46)  from  (4-30)  we  note  that  each  term  involving  w  as  a  negative 
exponent  -goes  to  zero  since  w  ->  oo  as  n  ->  oo .  Since  z0  is  a  constant  after 
many  passes,  the  first  term  in  (4-46),  which  represents  the  effect  of  the  initial 
optical  pulse  shape  on  the  pulse  shape  after  n-+<x>,  has  been  reduced  to  a 
constant.  This  indicates  that,  after  many  passes,  the  optical  pulse  shape  is  not 
affected  by  the  characteristics  of  the  initial  optical  pulse.  Therefore,  this  theory 
predicts  that  an  initially  random  optical  pulse  evolves  into  the  same  shape  as  an 
initially  smooth  optical  pulse.  It  follows  that  the  theory  should  be  applicable  to 
FELs  starting  from  classical  spontaneous  emission  through  coherence 
development  to  the  onset  of  strong  optical  fields.  It  should  be  noted  that  the 
spontaneous,  random  field  described  here  is  classical,  and  does  not  involve  the 
quantum  processes  analyzed  in  the  previous  section. 
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To  observe  the  evolution  of  coherence  development  in  a  more  quantitative 
way,  several  numerical  solutions  of  equation  (4-7)  were  obtained.  Shown  in 
Figure  4-5  is  the  coherence  development  of  an  initially  random  continuous  wave 
optical  field  with  average  strength  a0  =  1,  dimensionless  current  j  =  1.8,  initial 
phase  velocity  v",  and  parabolic  electron  pulse  shape  o(z)  =  1  -  2z2  /  az2  for 
\z\  <  J2az  and  zero  for  \z\>  V2az,  of  width  az  with  losses  determined  by  Q=6. 
The  desynchronism  is  cf  =  0.02,  and  the  initial  optical  amplitude  has  a  random 
component  added  at  each  size  z.  The  plot  of  y'(z-r)  at  the  lower-left  panel 
shows  the  parabolic  electron  pulse  shape  and  the  small  slippage  as  r  goes  from 
zero  to  one.  The  optical  amplitude  evolution  |a(z,r?)|  in  the  middle-left  panel 
shows  an  initially  random  pulse  which  becomes  much  more  coherent  over  n  =  40 
passes.  The  final  pulse  shape  in  the  upper-left  panel  is  much  smoother  and 
shows  a  high  degree  of  coherence.  The  gain  spectrum,  gr{v)  in  the  lower- 
center  panel  shows  the  antisymmetric  gain  function  for  reference.  The  evolution 
of  the  logarithm  of  the  optical  power  spectrum,  log  (P(v,n))  is  shown  in  the 
middle-center  panel.  The  evolution  of  coherence  leads  to  a  spectrum  centered  at 
the  phase  velocity  for  maximum  gain.  The  final  power  spectrum  in  the  upper- 
center  panel  is  narrow  indicating  long  range  phase  coherence.  The  optical 
phase  evolution  <f{z,n)  in  the  middle-right  panel  shows  the  development  of  the 
optical  phase  over  n  =  40  passes  with  the  final  phase  profile  at  the  upper-right. 
The  power  evolution  P(n)  in  the  lower-right  panel  shows  power  initially  decrease 
as  incoherent  parts  of  the  optical  field  are  removed  by  losses.  Then  the  power 
begins  to  grow  as  the  pulse  becomes  more  coherent  and  more  centered  at  the 
frequency  for  maximum  gain.  Figure  4-5  also  implies  that  the  optical  shape  is 
determined  primarily  by  the  shape  of  the  electron  pulse. 
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**  Analytic  FEL  Pulse  Evolution  ** 
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Figure  4-5.  Optical  coherence  development  from  a  random  pulse  solving  (4-7)  numerically. 


E.  CONCLUSIONS  AND  FURTHER  RESEARCH  SUGGESTIONS 

By  assuming  a  long  electron  pulse,  weak  optical  fields,  and  low  gain,  the 
electron  and  optical  pulse  shapes  have  been  expanded  in  a  Taylor  series  in 
order  to  solve  the  pendulum  and  wave  equations.  The  result  is  a  single  first- 
order,  partial-differential  equation  describing  the  evolution  of  an  optical  pulse. 
By  including  the  effects  of  desynchronism  and  resonator  losses,  a  more 
complete  partial-differential  describes  optical  evolution  over  many  passes  in  an 
FEL  oscillator.  As  a  direct  consequence  of  deriving  the  differential  equation,  the 
familiar  antisymmetric  gain  curve  and  symmetric  optical  phase  shift  are  obtained, 
as  well  as  functions  describing  the  effects  of  slippage  and  the  electron  pulse 
shape.    Solving  the  differential  equation  (4-7)  for  a  continuous  electron  beam 
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yields  the  simple  single-mode  gain  value  of  0.135;  and  results  in  no  change  in 
the  optical  pulse  shape.  Electron  pulses  with  parabolic  and  triangular  shapes 
are  shown  to  yield  gains  less  than  predicted  by  single-mode  gain  calculations 
since  there  are  fewer  electrons  distributed  over  the  entire  length  of  the  electron 
pulse.  The  shape  of  the  electron  pulse  is  shown  to  determine  of  the  shape  of 
the  optical  pulse  in  steady  state  after  many  passes.  Coherence  evolution  of  an 
initially  random  pulse  is  also  studied. 
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V.  SHORT  PULSE  EVOLUTION  IN  THE  BOEING 
APLE  OSCILLATOR 


A.  INTRODUCTION 

The  previous  chapters  have  looked  at  FELs  with  low  gain  and  weak  optical 
fields.  We  now  examine  an  FEL  oscillator  that  can  produce  high  gain  and  strong 
fields,  but  is  designed  to  produce  only  moderate  optical  fields  as  a  seed  laser 
for  an  FEL  amplifier.  Indeed,  the  problem  investigated  is  how  to  make  an 
oscillator  work  poorly,  or  to  reduce  the  electron-light  coupling,  so  that  only  a 
small  energy  spread  is  induced  in  the  electron  beam.  This  chapter  discusses 
the  Boeing/LANL  Average  Power  Laser  Experiment  (APLE)  which  is  designed  to 
produce  100  kW  of  average  power  for  3  minutes  from  a  FEL  using  a  Single 
Accelerator,  Master  Oscillator,  Power  Amplifier  (SAMOPA)  design,  shown  in 
Figure  5-1  [19,20]. 
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Figure  5-1 .  Schematic  diagram  of  the  Boeing  APLE  SAMOPA 
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In  Figure  5-1,  the  accelerator  section  (2)  preceding  the  FEL  oscillator  (3) 
gives  the  beam  17  MeV  energy  with  an  instantaneous  energy  spread  of  0.3%  at 
an  average  current  of  0.23  A  and  a  normalized  emittance  (90%  edge)  of 
e„  =  80 n  mm-mrad.  The  electron  beam,  with  a  peak  current  of  140  A  is  then 
directed  into  the  oscillator  (3)  where  it  produces  A  =  10.6/y/77  wavelength 
radiation  of  approximately  100  W  average  power.  Each  electron  micropulse 
may  have  an  energy  slew  of  up  to  3%  across  the  pulse.  Emerging  from  the 
oscillator,  the  electron  beam  is  accelerated  (4),  and  then  may  be  compressed  in 
two  interleaved  bunchers  (5)  to  increase  peak  current.  The  beam  is  then 
accelerated  further  (6)  to  34  MeV  and  the  energy  slew  removed.  The  electron 
beam  enters  the  power  amplifier  (7)  with  a  peak  current  of  450  A  and  an 
instantaneous  energy  spread  of  <  1.0%.  In  the  amplifier,  the  optical  pulse  from 
the  oscillator  experiences  a  gain  of  G  « 1000  to  meet  the  design  criteria  of  100 
kW  average  laser  power  for  3  minutes  [19,20]. 

This  chapter  examines  the  goal  of  passing  an  electron  beam  through  the 
oscillator  (3)  without  increasing  its  energy  spread  beyond  «1.0%.  This 
procedure  produces  an  optical  pulse  with  100  W  average  power,  and  determines 
an  optical  pulse  shape  and  size  that  optimizes  the  amplifier  efficiency.  The 
effect  of  an  energy  slew  in  the  electron  beam  is  investigated,  and  we  also 
consider  the  optimum  design  for  a  high  power  oscillator. 

B.  OSCILLATOR  DESCRIPTION 

The  APLE  oscillator  uses  a  linearly  polarized  undulator  of  wavelength 
A0  =  2.36  cm  over  N  =  100  periods  for  a  length  of  Z_  =  Ml0  =  2.36  m  and 
undulator  parameter  K  =  eBA0/27rmc2  =  0.23.  The  optical  wavelength 
produced       is       determined       from       the       resonance       condition       as 
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X  =  XQ\\  +  K2)J2y2  «  10.6/ym.  The  relatively  large  optical  mode  waist, 
W0  =  0.2cm,  compared  to  the  small  matched  electron  beam  radius  [1], 
re  =  yj(£n  ^o/2  nK)  =  008  cm,  leads  to  an  averaged  filling  factor  of 
F  =  (^e/^o)2/(1  +  l2/122o2)  m  0.12  with  Rayleigh  length  Z0  =  L/2.  The  electron 
beam  density  p=  1.4  x  1012  cm'3  gives  a  dimensionless  current  of 
7  =  8A/(  7ieKLJJ)2p/y3mc2  =250  where  the  Bessel  function  factor 
JJ  =  j0(£)  _  J^)  «  0.987  with  4  =  K2/2(l  +  K2)  «  0.025  is  present  in  order  to 
account  for  the  linearly  polarized  undulator  [1].  The  theoretical  single-mode  gain 
is  estimated  by  G  =  0.135y'F  a  58%,  but  is  not  self  consistent  and  does  not 
include  the  effects  of  energy  slew,  beam  quality,  or  optical  mode  distortion.  The 
length  of  the  electron  micropulse  is  /e  =  1.8  cm,  or  oz  «  17  slippage  distances, 
Nl. 

The  electron  phase  velocity,  v  =  L[{k  +  k0)pz  -/c],  measures  the  resonance 
between  the  electron,  undulator,  and  optical  fields  [1].  An  electron  beam  energy 
spread  causes  a  Gaussian  spread  in  z  velocities  and  phase  velocities  with 
standard  deviation  Av=  aG  =  AnN^yjy  »  1  where  Ay/y=  0.0015  for  the  APLE 
oscillator.  A  linear  energy  slew  of  3%  over  the  electron  micropulse  causes  a 
ramp  in  phase  velocities  from  v«  -25  to  v«  25.  An  angular  spread  of  the  beam 
causes  an  exponential  spread  in  the  phase  velocity  characterized  by  a  1/e  width 
of  Ave  =ae  =1.5.  The  energy  and  angular  spreads  do  not  significantly  degrade 
APLE  oscillator  performance.  The  large  energy  slew  of  Av«50  over  oz  «  17 
slippage  distances  means  that  each  slippage  length  of  the  beam  experiences  a 
phase  velocity  shift  of  Av«  3  during  each  pass. 
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C.  DIFFRACTION  EFFECTS 

When  the  effects  of  optical  diffraction  are  included  in  the  FEL  interaction, 
the  phase  velocity  for  maximum  gain  is  shifted  by  the  Rayleigh  length  z0  [1].  The 
electron  phase  velocity  that  gives  the  optimum  peak  single  pass  gain  was  found 
to  be  \%ax  w  2.6  +  1/z0  *  5  [1].  Optical  mode  distortion  results  from  high  current 
which  amplifies  the  radiation  overlapping  the  electron  beam.  This  effect  causes 
the  optical  mode  at  the  end  of  the  undulator  to  decrease  in  size  and  increase  the 
filling  factor.  A  single  pass,  three  dimensional,  self-consistent  numerical 
simulation  including  diffraction,  electron  beam  radius,  and  beam  quality  was 
used  to  determine  the  electron  phase  velocity  that  produces  the  peak  single- 
pass  gain,  vg73*  ,  and  the  new  effective  filling  factor.  The  simulation  used  a  peak 


current  density  of  j-250,  normalized  electron  beam  radius  <je  =  r^njLX  =  0.28 
with  a  parabolic  shape  [1],  energy  spread  with  standard  deviation  <jg  =1 
distributed  as  a  Gaussian,  angular  spread  described  by  an  exponential 
distribution  of  1/e  width  <je=1.5,  and  a  weak  optical  field. 

The  strong  fields  attained  over  many  passes  and  the  resonator  mirrors  can 
affect  the  steady-state  distortion  that  occurs  in  steady-state.  A  multiple-pass 
numerical  simulation  including  resonator  mirrors  and  losses  was  next  used  to 
determine  the  optical  mode  distortion  over  multiple  passes  in  strong  optical 
fields.  The  losses  were  characterized  by  Q=5,  where  1/Q  represents  the 
fractional  loss  of  optical  power  per  pass  in  the  resonator  from  both  outcoupling 
and  absorption.  The  resonator  was  characterized  by  a  mirror  radius  so  that  the 
fundamental  mode  has  edge  loss  around  the  mirrors  of  1%.  Strong  optical  fields 
cause  the  gain  to  decrease  until  the  gain  just  equals  the  loss  on  each  pass 
through  the  undulator.  Since  the  gain  is  reduced,  the  amount  of  mode  distortion 
is  reduced.    The  effective  filling  factor  now  obtained  by  comparing  the  cross 
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sections  of  the  undistorted  optical  mode  with  the  distorted  optical  mode  obtained 
from  the  multiple  pass  simulation  gives  F  «  0.04.  It  is  estimated  that  transverse 
diffraction  and  mode  distortion  in  the  FEL  oscillator  can  be  summarized  by  using 
an  effective  y'F«10  in  two  dimensional  codes  that  explore  longitudinal  pulse 
effects. 


D.  PULSE  EFFECTS 

Running  the  APLE  oscillator  with  a  relatively  large  current  of  jF  =  10,  and 
extracting  only  100W  average  power,  while  inducing  less  than  1%  energy 
spread  requires  the  electron-optical  coupling  to  be  significantly  reduced.  The 
first  two  cases  presented  below  examine  options  for  reducing  the  coupling.  The 
last  case  increases  the  coupling  to  achieve  high  efficiency  without  the  amplifier. 

1.  No  Energy  Slew  in  the  Undulator 

One  method  used  to  reduce  electron  -  optical  beam  coupling  is  to  operate 
at  a  large  desynchronism,  d  =  -2AS/NA,  where  AS  is  the  mirror  deviation  from 
exact  synchronism  [1].  Small  desynchronism  values  are  known  to  cause 
unstable  operation  because  a  small  change  in  resonator  length  can  alter  the 
FEL  from  high  power  to  low  power,  or  even  zero  power  [1].  Figure  5-2  shows 
optical  power  versus  desynchronism  for  the  APLE  oscillator.  It  was  obtained 
using  a  two  dimensional  numerical  simulation,  including  slippage  effects  with  an 
electron  pulse  length  of  crz«17,  a  Gaussian  spread  in  electron  phase 
velocities  of  oG  =  1.5,  and  outcoupling  and  losses  characterized  by  Q  =  5. 
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Figure  5-2.  Power  versus  desynchronism  for  the  APLE  oscillator. 


A  larger,  more  optimum  desynchronism  of  d  =  0.34  gives  low  power,  a  long 
optical  pulse,  and  a  small  induced  energy  spread,  as  shown  in  Figure  5-3.  The 
long  parabolic  electron  pulse  of  length  a2=M  is  displayed  in  the  lower-left 
panel  showing  the  slippage  over  one  pass  through  the  undulator  from 
dimensionless  time  r=  ct/L  =  0  to  x- 1.  The  middle-left  panel  shows  the  optical 
pulse  evolution  growing  from  shot  noise  over  n  =  800  passes  to  a  final  peak 
optical  amplitude  of  \a\  =  18  shown  at  the  top-left.  The  shot  noise  is  introduced 
by  including  a  random  electron  phase  with  standard  deviation  of  S£=  10"5.  The 
middle-center  panel  shows  the  evolution  of  the  optical  power  spectrum  over 
n  =  800  passes  with  the  top-middle  panel  displaying  the  final  power  spectrum. 
The  bottom-center  panel  shows  the  single-mode  gain  spectrum  for  reference. 
The  evolution  of  the  electron  distribution  over  the  n  =  800  passes  is  shown  in  the 
middle-right  panel  with  the  final  distribution  at  the  top-right.  The  bottom-right 
panel  displays  the  evolution  of  the  total  optical  power  within  the  oscillator  over 
/?  =  800  passes,  with  the  final  dimensionless  power  P  =  86  corresponding  to 
approximately  2200  W.     Since  the  trailing  half  of  the  electron  pulse  never 
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interacts  with  the  optical  pulse,  most  of  the  energy  distribution  is  kept  within  a 
narrow  band.  Only  the  electrons  on  the  leading  edge  of  the  electron  pulse  give 
up  energy  to  the  light  which  produces  the  small  bump  on  the  left  of  the  final 
electron  distribution  spectrum.  The  final  energy  spread  is  only  about  0.6%  and 
satisfies  the  APLE  criteria. 
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Figure  5-3.  Optimum  desynchronism  of  d=0.34  produces  required  power  without  inducing  large 

energy  spread. 


"The  effects  of  a  small  random  "jitter"  in  the  electron  energy  coming  from 
the  accelerator  was  also  explored.  The  study  used  the  two-dimensional 
simulation  but  included  a  small  random  phase  velocity  component  to  each 
electron  micropulse.  A  timing  jitter  was  included  by  adding  a  small  random 
change  to  the  desynchronism  length  on  each  pass.      These  jitters  had  a 


59 


significant  effect  on  small  values  of  desynchronism  and  caused  the  FEL  to  not 
operate  for  the  negative  desynchronism  values.  This  effect  could  be  expected 
since  the  region  of  small  desynchronism  in  Figure  5-2  is  considered  unstable. 
The  jitters  had  little  or  no  effect  on  the  power  at  moderate  to  large 
desynchronism  values. 

2.  Energy  Slew  in  the  Oscillator 

Another  option  for  the  APLE  SAMOPA  design  induces  a  linear  electron 
energy  slew  of  between  1%  and  3%  along  the  micropulse.  The  effect  is  to 
reduce  electron-optical  beam  coupling  in  an  attempt  to  leave  large  portions  of 
each  micropulse  unaffected  by  interactions  with  the  radiation  field.  The  energy 
slew  is  converted  to  a  slew  in  the  electron  phase  velocities  by  Av=  AnNkyjy. 
The  slew  in  phase  velocities  over  a  pulse  length  of  az  =  17  slippage  distances 
results  in  a  shift  of  n/4  (1%  slew)  or  3tt/A  (3%  slew)  over  each  slippage 
distance.  Superimposing  a  Gaussian  energy  spread  of  <rG  =  1  on  the  slewed 
pulse  causes  light  to  interact  with  electrons  with  phase  velocity  differences  of 
«;r.  Since  the  change  in  phase  velocity  is  greater  than  the  gain  bandwidth,  the 
FEL  interaction  is  affected.  In  order  to  reach  a  stable  power  and  stable  optical 
pulse  shape  when  there  is  a  significant  energy  slew,  the  desynchronism  must  be 
near  zero  so  that  the  light  continues  to  interact  with  the  same  energy  component 
of  each  electron  micropulse.  Unfortunately,  when  the  desynchronism  is  small, 
the  FEL  coupling  is  large  resulting  in  strong  optical  fields  and  a  large  electron 
energy  spread  at  the  end  of  the  interaction.  A  large  desynchronism  causes  the 
light  to  experience  a  significant  change  in  resonance  condition  as  discussed 
above.  In  that  case,  a  chaotic  optical  pulse  shape  and  power  evolution  are 
observed  in  simulations. 
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Figure  5-4  shows  the  results  of  a  two-dimensional  multimode  numerical 
simulation  with  jF  =  10,  a  large  desynchronism  of  d  =  0.06,  energy  spread  of 
erG  =1.5,  and  an  energy  slew  of  3%  over  n  =  1500  passes.  The  middle-left  panel 
displays  the  chaotic  optical  pulse  evolution  with  the  final  pulse  shape  shown  at 
the  top-left.  The  two  upper-center  panels  show  the  broad  optical  power 
spectrum  evolution  due  to  a  broad  range  of  resonance  conditions  contained  in 
the  slewed  electron  pulse.  The  lower-right  panel  displays  the  chaotic  power 
evolution  in  the  resonator  with  about  50%  power  variations  over  a  few  hundred 
passes. 
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Figure  5-4.  APLE  oscillator  with  3%  linear  electron  energy  slew  present. 

The  three  upper-right  panels  display  the  electron  energy  distribution  as  a  phase 
velocity  distribution  at  the  end  of  passes  n  =  1,  750,  and  1500.  On  the  first  pass, 
the  3%  energy  slew  and  the  0.3%  Gaussian  energy  spread  in  the  electron 
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micropulse  are  shown.  In  later  passes,  the  electron-light  interaction  induces 
two  to  three  times  the  energy  spread  originally  present,  which  is  the  maximum 
allowable  by  APLE  specifications. 

3.  High  Power  Oscillator  Experiment 

Another  possible  goal  for  the  APLE  laser  system  is  the  operation  of  the 
oscillator  alone  at  the  high  average  power  of  10  kW  for  3  minutes.  This  goal 
requires  a  complete  change  in  design  criteria  and  objectives  from  those  of  the 
previous  two  sections.  Rather  reducing  the  electron-optical  beam  coupling,  the 
goal  here  is  to  increase  coupling  to  get  maximum  power  from  the  FEL  oscillator. 

The  first  method  for  increasing  power  simply  increases  the  current  density 
j  ccK2  by  increasing  the  undulator  parameter  K.  This  is  accomplished  by 
reducing  the  gap,  g,  between  the  undulator  magnets,  since  K  oc  e'"9^0 .  The 
increase  in  coupling  is  expressed  by  y  occurs  because  increasing  Kwill  increase 
the  transverse  wiggling  motion  of  the  electrons;  the  amplitude  of  transverse 
oscillations  is  Ax  =  KA0/27ry. 

Secondly,  the  FEL  power  can  be  increased  by  using  a  small  value  of 
desynchronism  as  seen  in  Figure  5-2. 

The  two  methods  above  increase  FEL  coupling  and  should  lead  to 
saturation  in  strong  optical  fields.  Saturation  occurs  when  a  significant  number 
of  electrons  have  lost  enough  energy  to  move  through  the  gain  spectrum 
bandwidth  Ay/y»~y2N,  which  corresponds  to  a  phase  velocity  change  of 
Av=47zNAy/y*-2K.  Another  view  of  the  saturation  process  uses  the  FEL 
resonance  condition  A  =/l0(7  +  K2)l2y2 .  If  a  significant  number  of  electrons  lose 
enough  energy,  (y-"\)mc2,  the  optical  wavelength  no  longer  remains  in 
resonance  with  the  electrons. 
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The  resonance  condition  reveals  several  ways  to  go  beyond  the 
saturation  point.  As  electrons  lose  energy  along  the  undulator,  the  undulator's 
wavelength,  XQ,  or  the  parameter,  K,  could  be  decreased  so  that  resonance  is 
maintained  along  the  undulator.  An  undulator  designed  to  take  advantage  of 
this  concept  is  known  as  a  "tapered  undulator".  Conceptually,  reducing  X0  or  K 
provides  a  phase  acceleration  for  the  electrons  with  the  effect  of  simply  adding  a 
phase  velocity  to  the  pendulum  equation  (2-14)  describing  the  electrons' 
microscopic  motion  [1]: 

OO  0  £ 

£=  v-  S+\a\cos{£+  <f),  where  £=  \kQ(f)df  +  kz  +  cot ,      (5-1) 

o 

where   /c0(r)  =  2;r//l0(  r)   is  the  changing  undulator  wavenumber  along  the 
undulator.  Such  tapering  yields  a  phase  acceleration  of  [1] 

5=[_2dko(zl 
dz 

In  order  to  increase  coupling,  the  undulator  parameter  is  increased  to 
K  =  0.8.  The  three-dimensional  numerical  simulation  using  K=0.8  yields  a  new 
effective  filling  factor  resulting  in  jF  =  70.  A  small  value  of  desynchronism 
d  =0.001  enhances  electron-optical  beam  coupling  further.  The  phase 
acceleration  5  =  8x  continues  to  extract  energy  from  the  electrons  after 
saturation.  Figure  5-5  displays  the  results  of  a  two-dimensional  numerical 
simulation  using  losses  represented  by  Q  =  10  over  n  =  1000  passes.  The 
optical  pulse  is  modulated  by  the  trapped-particle  instability  and  grows  to  a  final 
peak  amplitude  of  \a\  =  570.  The  trapped-particle  instability  occurs  when 
electrons  trapped  in  phase  space  oscillate  at  the  synchrotron  frequency  leading 
to  the  formation  of  several  sidebands.    The  final  average  power  P  =  2.5x10" 
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corresponds  to  approximately  20  kW    and  is  twice  the  high  power  goal  of  the 
oscillator. 
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Figure  5-5.  High  average  power  from  tapered  oscillator. 

E.  CONCLUSIONS  AND  FURTHER  RESEARCH  SUGGESTIONS 

The  Boeing/LANL  Average  Power  Laser  Experiment  (APLE)  design  required 
limiting  the  electron-optical  coupling  in  the  oscillator  in  order  to  reduce  the 
energy  spread  of  the  electron  pulse  prior  to  injection  into  the  amplifier.  Several 
methods  of  limiting  the  coupling  were  examined,  including  the  use  of  large 
desynchronisms  and  the  use  of  an  energy  slew  in  the  electron  pulse.  Large 
desynchronisms  were  found  to  meet  the  requirements  of  a  smooth,  relatively 
long  optical  pulse  with  less  than  double  the  initial  energy  spread  of  the  electron 
pulse. 
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Using  an  energy  slew  to  reduce  electron-optical  coupling  resulted  in  a 
chaotic  optical  pulse  and  energy  extraction  from  a  large  portion  of  the  electron 
pulse. 

Producing  a  high  average  power  optical  pulse  using  a  tapered  oscillator 
without  the  amplifier  was  examined  in  this  chapter.  Further  investigation  on  the 
ability  of  the  mirrors  to  withstand  the  high  optical  powers  incident  on  them,  or  a 
proposal  for  a  mirror  system  that  can  operate  at  those  powers,  is  needed. 
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VI.  WHISTLER  WAVE  AMPLIFICATION  IN  THE 
MAGNETOSPHERE 


A.  INTRODUCTION 

1 .  Description  of  Whistler  Waves 

Whistler  waves  are  very  low  frequency  (VLF)  electromagnetic  radiation 
that  propagates  through  the  earth's  ionosphere  and  magnetosphere,  and  are 
detected  at  magnetically  conjugate  field  points  on  the  earth's  surface  [21]. 
Whistler  waves  were  first  commonly  detected  during  World  War  I  by  each  side 
as  they  used  amplifiers  on  long  telephone  lines  to  listen  to  the  other's 
conversations  [22],  Instead,  what  they  sometimes  heard  were  whistling  sounds 
decreasing  in  frequency  from  several  to  one  kHz  over  a  period  of  about  one 
second.  The  long  telephone  lines  were  acting  as  antennas  for  these  whistlers, 
and  the  high  gain  devices  boosted  the  otherwise  weak  signal.  It  also  found  that 
the  rate  of  whistler  activity  varied  with  time  of  day,  magnetic  activity,  and  location 
on  the  earth  [21].  Whistlers  are  more  common  at  night,  in  mid-latitudes,  and 
during  periods  of  increased  lightning  activity  [22]. 

The  correlation  between  periods  of  high  lightning  activity  and  whistler 
detection  led  to  the  following  explanation  for  their  occurrence.  A  lightning  flash 
in  the  southern  hemisphere  dumps  electromagnetic  energy  into  the  atmosphere 
in  a  short  broad-band  pulse.  This  pulse  may  then  propagate  through  the 
ionosphere  and  into  the  magnetosphere.  The  magnetosphere  is  highly 
dispersive,  so  as  the  pulse  travels  along  a  magnetic-field-aligned  duct  it  is 
spread  out  due  to  higher  frequencies  traveling  faster  than  lower  ones.  The  wave 
may  also  be  amplified  as  it  propagates  along  the  field  line  by  interactions  with 
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energetic  charged  particles.  As  the  wave  emerges  from  the  magnetosphere  at 
the  conjugate  point  in  the  northern  hemisphere,  it  is  detected  as  a  whistler 
beginning  with  higher  frequencies  and  trailing  off  to  lower  frequencies. 

Other  mechanisms  have  also  been  found  to  generate  whistler  waves. 
Morse  code  dashes  longer  than  150  milliseconds  may  produce  whistlers 
accidentally  [22].  VLF  transmitters  placed  in  mid-latitudes  may  stimulate 
whistlers  detected  at  the  conjugate  point.  Over  the  last  forty  years,  whistler 
waves  have  primarily  been  studied  using  VLF  transmitters  to  simulate  the 
injected  whistler  [22]. 

2.  Description  of  the  Magnetosphere 

The  immediate  environment  of  the  earth  consists  of  a  neutral  atmosphere 
extending  up  to  about  60  km  above  the  earth's  surface,  the  ionosphere,  a  region 
consisting  of  increasingly  ionized  gas  as  well  as  neutral  particles  extending  up  to 
about  1000  km  in  altitude,  and  the  magnetosphere,  consisting  of  ionized 
hydrogen  and  free  electrons,  extending  out  to  distances  of  about  100,000  km 
from  the  earth  [23].  There  is  no  distinct  boundary  between  the  ionosphere  and 
magnetosphere;  it  is  rather  a  smooth  continuum. 

The  structure  and  processes  of  the  magnetosphere  are  controlled  by  the 
earth's  magnetic  field  and  the  solar  winds  from  the  sun.  The  interaction  of  the 
solar  wind  and  the  earth's  magnetic  field  forms  the  bow-shock  at  approximately 
12Re  (where  Re  -  6370  km  is  the  mean  radius  of  the  earth)  with  the  solar  wind 
dragging  field  lines  around  the  earth  forming  a  tail  on  the  anti-solar  direction. 

The  region  in  space  in  which  the  earth's  magnetic  field  is  closely 
approximated  by  a  centered  dipole  field  inclined  at  11°  to  the  spin  axis  is  known 
as  the  inner  magnetosphere  [23].  The  inner  magnetosphere  extends  out  to 
about  7R  and  is  populated  by  a  background  'cold'  plasma  with  energies  of  0.1 
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to  1  eV  and  a  less  dense  'hot'  plasma  with  energies  from  1  keV  to  100  MeV  [23]. 
These  energetic  charged  particles  are  magnetically  trapped  in  helical  orbits 
spiraling  along  field  lines  and  make  up  the  radiation  belts.  Figure  6-1  shows  the 
inner  magnetosphere  with  dipole  magnetic  field  lines  and  the  helical  paths  of 
energetic  trapped  particles.  The  cold  plasma  is  responsible  for  the  index  of 
refraction  for  VLF  while  the  part  of  the  hot  plasma  that  is  resonant  with  the  VLF 
wave  can  amplify  the  wave.  A  typical  value  for  the  index  of  refraction  is  n  «20 
to  30  due  to  the  cold  plasma. 


Figure  6-1 .  Inner  magnetosphere  with  dipole  field  lines  and  energetic  particle  paths. 

Approximating  the  earth's  magnetic  field  as  a  centered  dipole  inclined  by 
11°  from  the  axis  of  rotation  yields  the  following  equation  for  the  magnetic  field: 
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B(R,A)  =  flJ^-j  J(l  +  sm2A),  (6-1) 


where  BQ  is  the  magnetic  field  at  the  surface  of  the  earth,  R  is  the  geocentric 
distance  to  the  field  line,  and  A  is  the  geomagnetic  latitude  [23].  A  dipole  field 
line  can  be  described  by 

R  =  Reqcos2A,  (6-2) 

where  R^  is  the  geocentric  distance  to  the  field  line  at  the  equator  [23].  Field 
lines  are  uniquely  identified  by  the  L  parameter 

L  =  -^-.  (6-3) 

R 

The  term  "ducting"  is  used  to  describe  the  bending  of  the  wave  as  it 
travels  up  along  a  field  line  in  the  southern  hemisphere  and  back  down  the  field 
line  in  the  northern  hemisphere.  Ducts  consisting  of  crests  of  higher  electron 
density,  or  troughs  of  lower  electron  density,  have  been  found  to  exist  along  field 
lines  in  the  magnetosphere  [22].  The  electron  density  gradient  creates  an  index 
of  refraction  gradient  within  a  duct  that  causes  a  wave  initially  entering  a  duct  to 
become  trapped  much  like  light  in  an  optical  fiber.  Since  the  wave  is  traveling 
from  a  higher  to  a  lower  index  of  refraction,  a  form  of  total  internal  reflection  may 
take  place  under  the  right  initial  conditions. 


B.  WAVE  PARTICLE  INTERACTIONS 
1.  The  Cold  Background  Plasma 

VLF  waves  propagating  through  the  magnetosphere  are  affected  by  the 
hot  and  cold  plasma  in  the  region  of  interaction.  The  electrons  in  the  cold 
plasma  are  several  orders  of  magnitude  less  energetic  than  the  electrons  in  the 
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hot  plasma  and  are  not  resonant  with  the  VLF  wave.  The  only  effect  of  the  cold 
plasma  on  the  VLF  wave  is  to  cause  dispersion.  The  cold  particle  density  pc  is 
several  orders  of  magnitude  greater  than  the  density  of  the  hot  particles. 

Since  the  VLF  wave  travels  along  a  magnetic  field  line,  a  useful 
coordinate  system  near  the  geomagnetic  equator  has  the  z-axis  follow  a  field 
line,  as  seen  in  Figure  6-2. 


vtl~ 


PARTICLE 


WAV* 


Figure  6-2.  A  field  line  showing  the  direction  of  propagation  of  the  wave  and  electrons  and  a 

relevant  coordinate  system. 


The  cold  particles'  motion  is  governed  by  the  Lorentz  force  law 

(  - 


mv  =  q 


ew+^(bw+b) 


(6-4) 


where  (..)=d(..)/dt,  q  is  the  charge  on  the  particle,  Ew  =EJ(  +  Eyy  is  the  VLF's 
arbitrary  electric  field,  Bw  is  the  VLF's  arbitrary  magnetic  field,  B=B0z  is  the 
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earth's  magnetic  field,  and  v  is  the  velocity  of  the  particle.  Since  the  free  ions 
are  much  more  massive  than  the  free  electrons,  they  can  be  considered  as 
being  immobile  at  VLF  driving  frequencies.  Assuming  the  wave  oscillates  at 
frequency  co,    the  linearized  equation  (6-4)  yields 


I 


-imcoVy  =  -e 


E,--£*. 


(6-5) 


and 


■imcovx  =  -e 


Ex  + 


Br 


(6-6) 


the  VLF  wave  magnetic  field  has  been  dropped  since  it  is  always  measured  to 
be  much  less  than  the  earth's  magnetic  field  [21  ].  Solving  equations  (6-5)  and 
(6-6)  algebraically  for  the  cold  electron  velocity,  we  find 


Vx  = 


moo 


/E, 


n 


CO 


( 


1- 


co2  J 


and 


(6-7) 


vy  = 


moo 


iE   +-EX 
V  oo 


co-  j 


-1 


(6-8) 


where  n  =  eB0/mc  is  the  electron  gyro-frequency. 

The  Lorentz  force  law  shows  how  the  VLF  wave  affects  the  motion  of  the 
cold  electrons.  The  wave  equation  is  required  to  see  how  the  cold  electrons 
affect  the  wave.  Ampere's  law  and  Faraday's  law  are 


vx4  =  £fjo+l^, 

c  c    at 


w         c    dt  ' 


(6-9) 
(6-10) 
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and  result  in  the  wave  equation 


\Anco  -,       co    =. 


-k(-k.Ew)  +  ^Ew=^Jc+^Ew,  (6-11) 


c2      c     c2    w' 
where  Jc  =  -epcv   is  the  cold  electron  current  density.     Substituting  the  cold 

electron  velocity  equations  (6-7)  and  (6-8)  into  the  wave  equation  (6-1 1 )  yields 

U2  -c2k2)Ey  =  — ^V-rfEy  +  AeA  (6-12) 


and 

(o>*-c>k2)Ex=       %   Jem-&E,\  (6-13) 

1  -  Q  ,  co   \  <o       J 

where  co2^  =  4ne2pc/m  is  the  cold  electron  plasma  frequency.  Solving 
equations  (6-12)  and  (6-13)  simultaneously  then  yields  the  dispersion  equation 
relating  k  and  co.  From  the  dispersion  equation  may  be  obtained  the  index  of 
refraction  caused  by  the  cold  electrons, 

co2  co(co±Q.) 

The  plus  sign  corresponds  to  the  index  of  refraction  for  a  left-hand  circularly- 
polarized  wave;  the  minus  sign  corresponds  to  a  right-hand  circularly-polarized 
wave.  The  whistler  wave  has  a  frequency  less  than  the  electron  gyro-frequency, 
and  must  be  a  right-hand  circularly-polarized  wave. 
2.  The  Whistler  Wave  Equation 

Having  developed  the  index  of  refraction,  which  is  interpreted  as  the 
effect  of  the  cold  electrons  on  the  whistler  wave,  the  effect  of  the  hot  electrons 
on  the  whistler  wave  is  now  explored.  Referring  back  to  Maxwell's  equations 
used  in  the  last  section,  the  current  density  must  now  include  the  hot  electrons 
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also:  J  =  Jh  +JC,  where  Jh  =  -evh±[-sin®,cos®,0]  is  the  single  particle  current 
density  and  0  =  J  Gl(t')dt'  is  the  phase  of  the  electron  cyclotron  motion.  When 
the  wave  envelope  is  a  narrow  band  whistler,  the  slowly  varying  amplitude  and 
phase  approximations  may  be  used  in  Maxwell's  equations,  dtEw«coEw, 
dzEw  «  kEw,  dtj> «  <o<t>,  dz<f> «  k<f>,  where  dz  =  d  {.)jdz  and  dt-d  {.)jdt  [1]. 
Having  shown  in  the  last  section  that  whistler  waves  are  right-hand 
circularly-polarized,  the  fields  can  be  expressed  as 

Ew  =  Ew  [cos  \f/,  -  sin  y/, 0] ,     Bw  =  Bw  [sin  y, cos  y/,  0] ,  (6-1 5) 

where  \\j  =  kz-ot  +$.     Inserting  the  whistler  fields  into  Ampere's  law  while 
allowing  the  amplitude  and  phase  to  evolve  results  in  the  following  equation: 

dzBw[-cos  y/,sin  y/,0]  +  Bw{k  +  dz(f)[sin  y/,cos  y/,0]  - 

1  E  4k-   (6"16) 

-dtEw[cos  y/,sin  y/,6]  +  — %-(-&+  dt(f>)[sin  y,cos  y/,0\  =  — J. 
c  c  c 

Defining  the  two  orthogonal  unit  vectors 

sy  =  [sin  if/, cos  y/,0]     and     e2  =  [cos  y/,-sin  y/,0]  (6-17) 

and  projecting  equation  (6-16)  onto  them  yields  the  following  two  scalar,  first- 
order  differential  equations: 

2     £ 

cBwdz<f>  +  Ew^  +  ckBw  -  coEw  =  -2-^-  -  4nevh±  cos(  y  +  0) ,       (6-1 8) 

£2-  go 

and 

cdzBw  +  dtEw  =  -4nevh±  sin{  y/  +  Q).  (6-1 9) 

These  equations  can  be  further  simplified  by  using  the  index  of  refraction 
equation  (6-14)  and  nEw  =BW  to  obtain 
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ncdzBz  +  dtBz  =  -AnenvhL  sin(  y/,  +  0) ,  (6-20) 

and 

Bw(ncdz(f>+dt4>)  =  -4nenvh±cos{y/+&).  (6-21) 

We  next  define  a  complex  whistler  field  b  as 

b  =  5§Le¥.  (6-22) 

mc 

Then  equations  (6-20)  and  (6-21 )  further  simplify  to 

'nc-?-  +  -§-)b  =  ^vexP(-itt+e)),  (6-23) 

\      oz     ox )  m 

where  v  =  nvh±/c  is  the  normalized  hot  electron  transverse  velocity  and 

£  =  kz  -  cri  +  \* Q.(t')dV 
is  the  longitudinal  electron  phase. 

The  single-particle  current  for  the  hot  electron  in  equation  (6-23)  is 
transformed  into  a  beam  of  electrons  by  a  weighted  average  over  their 
transverse  velocities  v,  their  transverse  phases  <9,  and  their  longitudinal  phase 
£  The  averaging  is  weighted  by  the  hot  electron  density  ph.  We  also  can 
assume  the  whistler  wave  is  spatially  flat  over  many  wavelengths  [26]  so  that 
d  bjd  z  «  0  to  obtain  the  whistler  wave  equation 

b  =  -ia>2ph(v-exP{-i(£+®))).  (6-24) 

A  whistler  wave  entering  the  magnetosphere  encounters  many  electrons 
randomly  spread  in  phase.  As  the  whistler  wave  equation  (6-24)  shows,  the 
interaction  with  randomly  spread  electrons  has  no  effect  on  the  whistler  and  the 
wave  does  not  grow.  However,  the  whistler  wave  affects  the  motion  of  the  hot 
electrons.  It  is  this  feedback  mechanism  that  allows  electron  bunching  and  wave 


74 


amplification.   The  effect  of  the  wave  on  the  motion  of  the  electrons  is  the  next 
topic  to  be  examined. 

3.  The  Hot  Electron  Equations  of  Motion 

The  motion  of  the  hot  electrons  in  the  presence  of  the  earth's  magnetic 
field  and  the  whistler  wave  fields  is  governed  by  the  Lorentz  force  law 


vh  =  - 


m 


ew  +  ^x{bw  +  b) 


(6-25) 
c 

However,  unlike  the  cold  electron  case,  the  whistler  magnetic  field  is  retained  to 

describe  electron  bunching.    Since  the  contributions  of  the  cold  electrons  has 

been  completely  contained  in  the  dispersion  relation  (6-14),  the  subscripts  "h" 

for  hot  and  "c"  for  cold  are  now  dropped  after  6-25. 

The  microscopic  bunching  of  the  electrons  due  to  the  wave-particle 

interaction  is  a  small  perturbation  on  the  macroscopic  gyro-motion  and  the 

bounce  motion  of  the  electrons  in  the  earth's  magnetic  field.    The  earth's  field 

near  the  geomagnetic  equator  is  not  exactly  uniform,  because  the  dipole's 

magnetic  field  lines  diverge  as  they  leave  the  north  pole  and  converge  as  they 

enter  the  south  pole.    This  results  in  a  component  of  the  field  in  the  radial 

direction  r ,  as  well  as  in  the  previously  defined  z  direction  [24].    In  order  to  find 

the  effect  of  the  radial  component  Br  of  the  earth's  magnetic  field  on  the  motion 

of  the  electrons,  Maxwell's  equation  V-B=0  is  used  to  find  6r: 


which  gives 


±-f(rB,)  +  ^  =  0  (6-26) 

r  dr  dz 


Br=iri^f]        •  (6"27> 

^  dz  Jr=0 
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The  Lorentz  force  in  the  z  direction  due  to  this  field  is 


1  ev±r  cBz 
c     dz 


mvz  =  -±^-^         ,  (6-28) 


where  v±  is  the  hot  electron's  gyro-velocity.    Averaging  over  one  gyration  and 
identifying  r  =  mv±c/eB  as  the  Lamor  radius  of  gyration  yields  [24] 


^  =  _^jS..  (6_29) 

z        2BZ  dz  v        ' 

Since  the  electron  energy  is  roughly  constant  while  spiraling  along  field  lines 

(vz2  +  v±2  =  constant),  we  have  vzvz  =  -v±v±,  so  that 

(6-30) 

Z.DZ     CZ 

and 


^x 

.  vxvz 

oBz 

2BZ 

dz  ' 

* 

^y 

V  V 

y  yv  z 

mz 

2BZ 

dz  ' 

(6-31) 

Recognizing  v]_/{2Bz)  in  equation  (6-29)  as  a  form  of  the  first  adiabatic  invariant 
[24],  it  is  therefore  a  constant.  Thus  equation  (6-29)  can  be  written  in  terms  of 
the  values  of  the  initial  transverse  velocity  v±0  and  magnetic  field,  BQ,  at  the 
equator. 

Adding  equations  (6-29),  (6-30)  and  (6-31)  into  (6-25)  gives  the  hot 
electron's  equation  of  motion  including  the  effects  of  the  whistler  wave  and  the 
earth's  inhomogeneous  magnetic  field.  Expressing  each  velocity  component 
yields 


e 
m 


Ex+UvyBz-vzBy) 


v^v7  SB 
+    x  z°       z ,  (6-32) 

2BZ    dz  v 
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vy  =  ~ 


m 


Ey+UvzBx-vxBz) 
c 


vyvz»  dBz 
2B7    dz  ' 


(6-33) 


Vz  =  ~ 


mc 


(VxB       vBx)-±3-**. 

v  x  y       y  xJ      2B0    dz 


(6-34) 


Since  the  region  of  space  most  often  identified  with  the  whistler  wave- 
particle  interaction  is  within  a  few  degrees  of  the  geomagnetic  equator,  a  useful 
approximation  to  the  earth's  dipole  magnetic  field  in  that  region  is 


Bz  =  B(z)  =  flb(l  -  qz2) , 


(6-35) 


where  q  =  9/(2L2Re2),  B0  is  the  magnitude  of  the  earth's  magnetic  field  on  the 
field  line  at  the  equator,  and  z  is  distance  measured  along  a  field  line  from  the 
equator  [25]. 

Substituting  the  form  in  equation  (6-35),  for  Bz  as  well  as  for  Bx  and  By 
and  the  right-handed  circularly-polarized  whistler  wave  from  equation  (6-15)  into 
equations  (6-32),  (6-33),  and  (6-34)  yields 


vx 

= 

_e/ 
m{ 

Vy 

= 

vz 

= 

_e_r 
~m{ 

Ecos  y/  +  ^-B(z)  — LnEcos  y/   +  vxvZoqz, 

CO  J 

(  v     i  \     v  A 

-Esin  \f/ — -B{z)  — -nEsin  y/   +  vyv^qz, 
c  c  J 

-2-nEcosy/ — -  nEsin  y/\-  (v±o)qz. 


(6-36) 
(6-37) 
(6-38) 


These  equations  can  be  further  simplified  in  terms  of  the  perpendicular  velocity 
vu  and  the  angle  0  between  the  whistler  wave  fields  and  the  electron's  velocity, 
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vy  =  v±  cos  0 ,  vx  =  -  v±  sin  0 .  (6-39) 

where  0  =  6+  J  Cl(t')dt' .   After  some  algebra,  the  equations  of  motion  reduce 
to 


(z-  -Vj^qz  +  -    "J-  sm( y  +  Q) , 

sin(  ^  +  0), 


mc 

eBw  ( .     nvz 
0  r?m  V        c 


.      o8(z)  .    eBM 


0  = 


+ 


mc       nmv± 


1- 


nv, 


\cos(y/+&). 


(6-40) 
(6-41) 

(6-42) 


Introducing  the  whistler  wave  amplitude  \b\  =  eBw/mc  and  the  electron  gyro- 
frequency  Q  =  eB(z)/mc  yields 


vz  =  -v]_  qz  +  \b\v±  sin{  y/  +  0) 


v±  =  vzv±Qz  + 


cb 


n 


1- 


nv, 


0  =  D(z)  + 


ebb 


A7V, 


1- 


A?V, 


S/f7(^  +  0), 


cos(y/  +  0) 


(6-43) 
(6-44) 

(6-45) 


The  first  terms  on  the  right  side  of  each  equation  (6-43)  through  (6-45)  describe 
the  fast  macroscopic  electron  motion  in  the  large  scale  dipole  field  of  the  earth. 
The  last  terms  are  proportional  to  the  VLF  wave  field  \b\  and  describe  the  small 
microscopic  scale  motion  in  the  combined  dipole  and  VLF  fields. 

Since  the  hot  electrons  spend  much  less  time  in  the  interaction  region 
than  the  bounce  period  in  the  dipole  field,  the  electron  motion  in  the  first  term  on 
the  right  side  of  equations  (6-43)  through  (6-45)  can  be  described  by 


■macro 


(t)  -  vzt. 


(6-46) 
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According  to  experimental  observations,  most  whistler  activity  and  VLF 
wave  amplification  occurs  when  co  =0.5/2,  known  as  the  "cutoff  frequency"  [22]. 
At  the  cutoff  frequency,  the  phase  and  group  velocities  of  the  wave  are  given  by 
vp  =  c/n  and  vg=dco/dk  =  c/n  respectively.  Since  at  resonance  v«0,  the 
electrons  satisfy  kv2o  =  co-Q  so  that  v2o  =  -  co/k  =  -vp  =  -  c'/n.  Hence,  the  wave 
at  the  cutoff  frequency  and  the  resonant  electrons  are  traveling  with  equal 
velocities  in  opposite  directions. 

Recognizing  that  0  has  both  fast  and  a  slowly  rotating  terms  allows  us  to 
separate  out  the  simple,  fast  periodic  gyro-motion  of  the  electrons  in  order  to 
concentrate  on  the  more  interesting  slower  evolution  of  the  bunching  electrons. 
Redefining  O  =  0{  +6S,  with  6f  =  fn(f ')<#',  allows  the  use  of  0s(f)  for  the  slow 

JO 

evolution  of  the  VLF  wave. 

The  arguments  within  the  trigonometric  functions  in  equations  (6-43) 
through  (6-45)  can  be  rewritten  in  terms  of  the  longitudinal  electron  phase 

£=  kz  -  cot  +  j'n(t')dt' ,  (6-47) 

where  ^  +  0  =  y/+  0f  +  0S  =  £+  0S  +  <f>.  Defining  the  electron  longitudinal  phase 
velocity  as 


v  =  C,  -  kvz  -  co  +  QQ 


<      qc2t2^ 

1  +  2—2- 

n 


(6-48) 


where  (6-46)  is  used  to  replace  z  in  terms  of  t  in  the  inhomogeneous  term.  Thus, 
the  electron  phase  acceleration  becomes 


3 

1=  y=kvz+4q^Yt.  (6-49) 
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Substituting  the  electron  phase  and  electron  phase  velocity  into  equations 
(6-43)  through  (6-45)  while  using  (6-46)  for  the  macroscopic  z  terms,  vZo  =  co/k 
for  the  macroscopic  terms,  and  v2o  =  -c/n  in  the  microscopic  terms  yields 


3 

v=  (a  -  vl)^-t  +  \b\va>sin{£+  0+  (f>),  (6-50) 


v  =  y^-t  +  2\b\sin(£+  0+  j),  (6-51) 

k 


0=  J-Acos(£+  o+  ^  (6.52) 

v 

where  the  definition  v  =  kvx  is  a  normalized  perpendicular  velocity,  v0  =  kv  is 
the  normalized  perpendicular  velocity  at  the  equator,  and  the  subscript  s  is 
dropped  from  #s,  since  all  coordinates  are  now  slowly  evolving.  The 
macroscopic  z  motion  along  the  inhomogeneous  field  line  uses  the  substitution 
z-ct/n.  The  microscopic  z  describing  electron  bunching  in  the  VLF  wave  uses 
the  fact  that  the  wave  and  electrons  are  counter-propagating  which  causes 
vz0  =  -c/n. 

The  first  term  on  the  right  side  of  equations  (6-50)  and  (6-51)  represents 
the  acceleration  of  the  hot  electrons  due  to  the  inhomogeneity  of  the  earth's 
magnetic  field.  The  acceleration  term  is  analogous  to  the  "taper"  term  in  a  free 
electron  laser  as  discussed  in  Chapter  V,  and  so  the  acceleration  term  is 
referred  to  as  a  "taper"  here.  Define 

S  =  — =-  =  — 5-5 — 5-,  (6-53) 

k2       2k2L2R* 

where  L  defines  the  number  of  geocentric  earth  radii  at  the  equator  to  the 

magnetic  field  line  being  used,  and  Re  =6370  km  is  the  mean  radius  of  the 
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earth.    Presenting  the  hot  electron  equations  of  motion  using  the  taper  term 
gives 

v=  (4  -  vl)Scot  +  \b\vcosin(£+  6+  <p)  (6-54) 

v  =  vSt  +  2\b\sin(£+  6+  <f)  (6-55) 

0=  -rlcos{£+  0+<f).  (6-56) 

v 

Several  interesting  observations  can  be  made  about  the  equations  of 
motion  (6-54)  through  (6-56).  Equation  (6-54)  follows  the  longitudinal  evolution 
of  the  electrons,  and  includes  the  pendulum  equation  with  a  phase  acceleration 
identified  as  taper.  The  taper  term  in  equation  (6-54)  shows  that  certain 
electrons  with  an  equatorial  normalized  perpendicular  velocity  of  v0±2 
experience  no  effect  on  their  longitudinal  phase  velocity  due  to  the 
inhomogeneity  of  the  earth's  magnetic  field.  A  normalized  perpendicular  velocity 
of  v0  =  ±2  corresponds  to  a  pitch  angle  of  «  ±63°.  The  electron  phase  velocity, 
v=  kvz{t)-  a>  +  Q(t),  shows  that  the  parallel  velocity  and  the  gyro-frequency 
are  each  changing  in  opposite  directions  at  the  same  rate  at  this  particular  pitch 
angle,  thereby  canceling.  The  perpendicular  velocity  is  only  slightly  affected  by 
taper  since  (6-55)  evolves  much  slower  than  the  longitudinal  phase  velocity 
equation  (6-46)  due  to  the  presence  of  the  large  factor  co  «  1 04/s. 


C.  NUMERICAL  SIMULATIONS 

The  equations  of  motion  (6-54)  through  (6-56),  together  with  the  wave 
equation  (6-24),  are  now  numerically  integrated.  The  solutions  reveal  many  of 
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the  same  characteristics  found  in  free  electron  lasers,  including  the  observation 
of  electron-bunching,  saturation,  trapped  particle  instability,  and  taper. 

The  parameters  used  for  the  simulations  are  taken  from  measurements 
made  by  the  Space,  Telecommunications,  and  Radioscience  Laboratory  at 
Stanford  University  from  Siple  Station,  Antarctica.  This  data  was  obtained  using 
their  >100  kW  transmitter  to  send  VLF  signals  at  between  one  and  six  kilohertz 
using  a  long  horizontal  dipole  antenna  [26].  The  amplified  signals  were  then 
received  at  the  conjugate  point  at  Lake  Mistissini,  Quebec  [26].  The  data  being 
simulated  here  is  a  transmitted  signal  at  co  =  104/s  estimated  by  the  observers 
through  time  of  flight  measurements  as  traveling  along  approximately  the  L  =  5 
field  line.  An  index  of  refraction  within  the  interaction  region  close  to  the  equator 
was  estimated  to  be  approximately  n  =  25.  The  hot  electrons  present  are 
represented  by  electrons  uniformly  spread  in  the  electron  phase,  £0  =  0  -»  2x, 
and  the  electron  cyclotron  phase,  60  =  0  -+  In.  The  time  of  interaction  is  taken 
to  be  7=1  s. 

The  first  case  examined  is  in  the  weak-field,  high-gain  regime.  In  an  FEL,  a 
weak  optical  field  is  identified  by  \a\<n.  Similarly,  a  weak  whistler  field  is 
identified  when  \b\coT2  <  n  which  corresponds  to  \b\  «  3  x  lO^/s  [27].  Figure  6-3 
shows  the  results  of  high  gain  in  a  weak  whistler  field  of  initial  value 
b0  =  \b(0)\  =  5  x  10"5/s  as  it  interacts  with  a  hot  plasma  of  frequency  of 
cob  =  0.2/s.  The  electron  beam  has  a  small  Gaussian  spread  in  phase  velocities 
characterized  by  width  standard  deviation  ov  =  0.5  /s  about  the  initial  phase 
velocity  v0  =0.  All  electrons  have  a  normalized  perpendicular  velocity  of  v0  =  1 
corresponding  to  an  initial  pitch  angle  of  x/4.  The  whistler's  frequency  is 
co  =  104  /s  with  an  interaction  time  of  7=1  second.  The  taper  term  is  S  =  0, 
which  indicates  a  homogeneous  magnetic  field.    An  FEL,  with  current  density 
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j>n  is  in  the  high  gain  regime  and  experiences  exponential  gain.  Similarly, 
when  o)2hcoT3  >  n,  (or  coh  >  -\nj(oTz  «  0.02/s)  the  whistler  wave  is  said  to  be  in 
the  high  gain  regime  where  exponential  growth  occurs  [27].  After  a  characteristic 
time  elapses  allowing  the  electrons  to  bunch  (approximately  0.35s  in  this 
example)  the  gain  becomes  exponential.  In  an  FEL  the  exponential  gain  is 
calculated  from  FEL{r)  «  exp|(y'/2)1'3V3r|/9  [1].  Similarly,  the  exponential 
gain  for  the  whistler  wave  is  calculated  by  [27] 


G{T)  *  exp 


Xatfillf&T 


/  9.  (6-57) 


Equation  (6-57)  yields  a  final  gain  after  the  one  second  interaction  of  G  ^34  dB. 
This  value  is  slightly  larger  than  that  shown  in  Figure  6-3  since  equation  (6-57) 
does  not  take  into  account  gain  degradation  due  to  a  spread  in  electron 
longitudinal  phase  velocities  or  the  delay  due  to  the  bunching  time.  The  phase 
is  shifted  linearly  by  nearly  the  same  factor  used  in  the  exponential  to  calculate 
the  gain,  A<p  «  (0^/2)  7/2  «  3  [27].  The  phase  space  in  (£  v)  displays  the 
final  positions  of  2000  sample  electrons  after  the  one-second  interaction.  Those 
electrons  going  up  in  phase  velocity  remove  energy  from  the  wave,  while  the 
electrons  going  down  in  phase  velocity  transfer  more  energy  to  the  wave  and 
cause  amplification.  The  bunching  of  electrons  in  phase  is  clearly  present  at  two 
locations  in  phase  space  since  there  are  two  phase  factors  (£,  =  0-»2;r  and 
00  =  0-»  2;r)  within  the  trigonometric  functions  in  the  equations  of  motion.  The 
two  dark  horizontal  lines  on  the  phase-space  plot  show  the  maximum  height  of 


the  closed-orbit  region  given  by  Avs  *  AJ\b\coT2  (see  Figure  6-3) 
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Figure  6-3.  Whistler  interaction  in  a  weak  VLF  field  with  high  gain. 

Most  whistlers  reach  strong  field  saturation  during  the  interaction  time. 
Once  saturation  occurs  the  electrons  become  trapped  in  closed  orbits  in  phase 
space  and  execute  synchrotron  oscillations  at  a  frequency  cos  *  ^jco\b\  [27]. 
These  oscillations  can  evidence  themselves  as  oscillations  in  the  power  and 
gain  of  the  whistler  wave.  Figure  6-4  shows  the  results  of  a  Gaussian 
distribution  of  hot  electrons  of  width  0^  =  10  /s  about  the  resonant  phase 
velocity  v0  =0  interacting  with  a  strong  whistler  field  for  a  time  of  T  =  0.4  s.  The 
electrons  are  at  a  pitch  angle  of  x/4  determined  by  v0=1.0,  with  a  density 
corresponding  to  a  hot  plasma  frequency  of  a>h  =  10/s.  The  whistler  wave  has 
an  initial  amplitude  b0  =  |b(0)|=0.01/s  with  frequency  co  =  104  /s  .  Again,  there  is 
no  taper  so  5  =  0.  The  phase  space  plot  of  Figure  6-4  shows  the  final  positions 
of  3000  sample  electrons  showing  their  bunching  at  two  phase  regions  and  each 
with  lower  energy.    The  spread  in  phase  velocities  causes  the  bunching  to  be 
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more  diffuse  and  creates  much  of  the  random  scatter  of  electrons  seen  in  the 
phase-space  plot.  Electrons  appearing  outside  the  final  separatrix  lines  have 
been  trapped  by  the  previously  larger  whistler  field  in  the  oscillating  pattern  seen 
in  the  gain  plot.  Note  that  the  final  gain  is  at  a  minimum  caused  by  the 
synchrotron  oscillations. 
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Figure  6-4.  Whistler  interaction  in  a  strong  VLF  field  with  high  gain. 


The  beginning  of  the  gain  curve  shows  that  after  the  short  bunching  time,  the 
gain  is  exponential  until  strong  field  saturation  occurs  at  a  value  of 
\b\  «  2(a>4h/Aco)  m  1.3  /s.  The  strong  fields  trap  the  electrons  and  cause  them 
to  perform  synchrotron  oscillations.  The  presence  of  approximately  4.5 
synchrotron  oscillations  is  clear  in  the  gain  plot  and  agrees  closely  with  the 


calculated  value  of  cos  «  yjco\b\  a  36 ;r  /s  which  would  predict  7  full  oscillations 
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over  the  interaction  time  of  0.4s.  The  small  discrepancy  between  the  number  of 
synchrotron  oscillations  observed  and  the  estimate  cos  «  36 ;r  /s  occurs  because 
of  the  spread  in  electron  phase  velocities.  Saturation  in  strong  fields,  electron 
trapping  and  synchrotron  oscillations  are  all  nonlinear  effects  observed  in  both 
FELs  and  whistler  waves  [27]. 

The  previous  figures  used  a  short  time  of  interaction  and  a  homogeneous 
magnetic  field  with  £=0.  For  a  longer  interaction  time,  the  effects  of  the 
inhomogeneity  in  the  earth's  magnetic  field  acts  as  a  phase  acceleration 
analogous  to  taper  in  an  FEL.  In  Figure  6-5,  the  electrons  begin  away  from  the 
geomagnetic  equator  and  far  off  resonance  with  a  Gaussian  spread  of  phase 
velocities  characterized  by  a  width  gv  =  1 0  /s  centered  at  vQ  =  2360  /s.  The 
electrons  have  an  equatorial  pitch  angle  of  njA  determined  by  v0  =  1,  with  a 
density  creating  a  hot  plasma  frequency  of  coh  =  1 0  /s.  These  hot  electrons 
interacting  with  an  initially  weak  whistler  field  of  b0  =  \b(0)\  =  0.01  /s  of  frequency 
Co  =  104  /s  for  a  time  of  T  =  1  s.  The  taper  value  of  8=  Q.2n  corresponding  to 
the  duct  at  L«5  with  index  of  refraction  n«25.  The  phase-space  diagram 
displays  the  final  positions  of  3000  sample  electrons.  The  electrons  begin  far 
above  resonance,  but  their  phase  velocity  decreases  as  they  move  into 
resonance  near  the  geomagnetic  equator.  The  simulation  is  designed  so  the 
sampled  electrons  are  in  resonance  with  the  wave  at  the  geomagnetic  equator  at 
a  time  f  =  7/2.  After  strong  coupling  between  the  hot  electrons  and  the  VLF 
wave  near  the  equator,  the  taper  causes  the  electrons  to  increase  their  phase 
velocity  as  they  move  past  the  equator  and  away  from  resonance.  However,  as 
seen  in  Figure  6-5  most  of  the  electrons  are  trapped  by  the  strong  whistler  wave 
field  and  continue  to  give  energy  to  the  wave. 
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Figure  6-5.  Whistler  interaction  in  a  inhomogeneous  magnetic  field  with  a  strong  VLF  field  and 

high  gain. 


The  gain  evolution  in  this  figure  shows  no  coupling  between  the  electrons  and 
whistler  until  the  electrons'  phase  velocity  is  nearly  resonant  at  time  t  =  0.5  s. 
Once  sufficiently  near  resonance,  the  electrons  become  bunched  and  the  gain 
becomes  exponential.  Synchrotron  oscillations  begin.  The  effect  of  the  taper 
causes  the  trapped  electrons  to  continue  to  give  up  energy  to  the  wave  after  the 
normal  saturation.  Most  of  the  electrons  remain  trapped  near  resonance  by  the 
strong  whistler  field,  even  though  the  earth's  dipole  field  provides  a  phase 
acceleration  to  move  them  away  from  resonance.  Those  electrons  not  trapped 
by  the  strong  VLF  wave  are  moved  off  resonance  by  the  taper.  The  gain 
evolution  shows  that  because  the  electrons  are  far  off  resonance  initially,  there 
is  no  growth  in  the  wave.  As  the  electrons  approach  resonance  at  the  equator, 
the  same  exponential  growth  seen  in  the  homogeneous  case  occurs  in  the  wave, 
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leading  to  saturation  and  synchrotron  oscillations.  The  inhomogeneity  of  the 
earth's  magnetic  field  provides  a  phase  acceleration  that  allows  continued 
energy  extraction  from  the  trapped  electrons.  This  is  seen  as  the  gain  slowly 
increases  even  after  the  homogeneous  saturation  level.  The  wave's  phase 
initially  decreases,  a  characteristic  of  tapering,  and  then  begins  to  grow  quickly. 
The  rate  of  phase  growth  decreases  in  conjunction  with  decreasing  gain  growth. 
The  general  characteristics  of  the  gain  and  phase  evolutions  shown  in  Figure  6- 
5  are  in  agreement  with  those  obtained  by  Carlson  et.  al.  of  Stanford  University 
[28].  Their  simulations  also  show  gain  initially  growing  exponentially  and  then 
reaching  either  a  steady  value  or  growing  slowly  [28].  They  also  found  the  phase 
would  initially  decrease  during  the  exponential  growth,  and  then  would  increase 
parabolically  [28]. 


D.  CONCLUSIONS  AND  FURTHER  RESEARCH  SUGGESTIONS 

Although  apparently  an  unrelated  process,  the  amplification  of  whistler 
waves  is  shown  in  this  chapter  to  be  closely  related  to  the  amplification  of  an 
optical  wave  in  a  free  electron  laser.  The  radiation  field  in  both  whistlers  and 
FELs  couples  in  resonance  to  "wiggling"  electrons  for  efficient  energy  transfer. 
Nonlinear  effects  such  as  saturation  and  the  trapped-particle  instability  are 
observed  in  both.  For  the  whistler,  the  inhomogeneity  of  the  earth's  magnetic 
field  lines  acts  in  an  analogous  manner  to  tapering  an  undulator.  The  electron 
equations  of  motion  and  the  wave  equation  derived  in  the  chapter  are  used  show 
many  of  the  same  effects  seen  in  whistler  and  FEL  research.  The  results  closely 
reproduce  the  results  of  other  researchers  using  more  traditional  approaches  to 
the  study  of  whistler  wave  amplification. 
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This  chapter  solves  for  whistler  wave  amplification  of  a  single  mode  using 
either  monoenergetic  hot  electrons  or  electrons  in  a  Gaussian  distribution  of 
phase  velocities  with  a  single  pitch  angle.  Further,  all  electrons  are  forced  to  be 
in  resonance  at  the  equator,  or  at  time  f  =  7/2,  regardless  of  their  initial  phase 
velocity.  The  research  needs  to  be  extended  by  adding  a  more  realistic 
distribution  of  hot  electron  energies  and  pitch  angles.  A  multimode  numerical 
simulation  should  be  developed,  similar  to  those  used  in  FEL  research,  so  that 
the  experimentally  observed  sweeps  in  frequency  of  the  whistler  waves  can  be 
investigated  numerically. 
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